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We derive a non-equilibrium fluctuation-dissipation relation for bosonic correlation functions from 
holography in the classical gravity approximation at strong coupling. This generalizes the familiar 
thermal fluctuation-dissipation relation in absence of external sources. This also holds universally 
for any non-equilibrium state which can be obtained from a stable thermal equilibrium state in 
perturbative derivative (hydrodynamic) and amplitude (non-hydrodynamic) expansions. Therefore, 
this can provide a strong experimental test for the applicability of the holographic framework. We 
discuss how it can be tested in heavy ion collisions. We also make a conjecture regarding multi-point 
, holographic non-equilibrium Green's functions. 
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I. INTRODUCTION 



Nature challenges us to understand phenomena in real time. This is very difficult to do in the microscopic formu- 
lation of the laws of nature given in the framework of quantum field theory. 

Though quantum field theory is experimentally successful in understanding microscopic processes when perturbation 
theory works, it has very limited success in describing macroscopic phenomena in real time. The best we can typically 
do is to calculate rate and cross-sections using perturbation theory, and use the results as inputs for kinetic or 
phcnomenological equations. As for example, for dilute quantum gases we can use the Boltzmann equation where 
the two body scattering cross-section is an input. However, for few systems we know how to formulate systematic 
corrections to the Boltzmann equations, which takes into account the uncertainty principles of quantum dynamics. 

The main difficulty is that time-dependent perturbation theory gives us behavior in time only in a Taylor series. 
This fails to give a uniform approximation in time away from the initial phase of evolution, and describe decoherence, 
thcrmalization and hydrodynamics in a unified framework. Even when the coupling constants are small, we need an 
essentially non-perturbative approach to understand the origin of irreversibility from the microscopic quantum field 
theory. 

In recent decades, there has been impressive progress in non-equilibrium quantum field theory, which attempts 
to bridge this gap using the two-particle irreducible action formalism 1 . Still a lot of progress is needed to develop 
approximation schemes which takes into account unitarity and conservation laws in a controlled manner. As these 
methods are non-perturbative, some progress is also needed in understanding renormalizability in these formalisms. 

Holography in the form of gauge/gravity duality is a new tool which gives a non-perturbative reformulation of 
quantum field theory Qj. This is particularly tractable when the quantum field theory is a strongly coupled gauge 
theory and the rank of the gauge group is large. In such a case, gauge/gravity duality maps a quantum field theory to 
a familar classical theory of gravity in one higher dimension. At present this is the best non-perturbative reformulation 
of quantum gauge theories at hand, which can describe phenomena in real time. 

In this paper, we will try to unearth a special feature of holographic duality at strong coupling and large rank of 
the gauge group (i.e. large N), which gives a genuinely non-equilibrium result and is also a model- independent (i.e. 
universal) feature 2 . By a genuinely non-equilibrium result, we imply one which cannot be obtained readily from 
equilibrium correlation functions. We also try to focus on a feature which gives fundamental insight on thermalization 
and nature of quantum kinetics. This leads us to study non-equilibrium fluctuation-dissipation relation as described 
below. 

It is known from non-equilibrium quantum field theory, that in order to get quantum formulation of kinetic equations 
we need the non-equilibrium spectral and statistical functions. We define them below. 

The spectral function can be understood as an off-shell generalization of density of states. Here we will define it 
for a bosonic operator 0(x,t). It can be obtained as a Wigner transform of the commutator, i.e. Fourier transform 
of the ralative coordinate as below : 



A(u),k,x,t) = J 
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The expectation value value above is taken in a non-equilibrium state. In case of fermionic fields we need to take the 
anti-commutator above. 

It can be proved (see appendix A) that the spectral function is proportional to the imaginary part of the Wigner 
transformed retarded Green's function as below : 



A(u),]&, x, t) = -2ImGjj(u,k, x, i). 



(2) 



Therefore the spectral function is real. In case of fermionic Hcrmitean operators it is also non-negative. 

Physically Aiw, k, x, t) encodes the probability that a quasiparticle with momentum k at a given point in space-time 
given by x and t will have energy u>. We note that the uncertainty principle forces us to consider the quasi-particle 
off-shell. At strong coupling we expect broadening of peaks and typically no sharply defined quasi-particles. 



1 This started with the closed-time-path formalism of Schwinger and Keldysh pF. The next advancement came with the two-particle 
irreducible (2PI) action formalisrn due to Cornwall, Jackiw and Tomboulis |2|; and incorporation of the Schwinger-Keldysh closed time 
contour in it by Calzetta and Hu H . This will be briefly reviewed in section flof this paper. For a review of recent progress of systematic 
approximations of the two-particle irreducible action and it's applications please see l4-c[. 

2 At this stage we can make a comarison to the famous universal holographic result at strong coupling and large TV that the viscosity to 
entropy density ratio r)js equals 1/4-7T fel. However, this is not a genuinely non-equilibrium result by our definition, because it can be 
computed via holographic prescritions from thermal correlators, and applying Kubo formula H. 
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At equiibrium, due to translational invariance, the spectral function A(ui, k, x, t) do not depend on the center-of- 
mass coordinates x and t. It is measurable by angle-resolved photoemission spectroscopy (ARPES) of electrons. Away 
from equilibrium, time resolved version of the same (tr- ARPES) is necessary (see for instance |10|). 

The statistical function can be understood as an off-shell generalization of quasi-particle distribution function in 
phase-space. It is the Wigner transform of the anti-commutator for bosonic fields (commutator for fermionic fields) 
as below : 

Gx:(w,k,x,t) = -i J d 3 r^e^- kr )({0(x+|,t+|),0(x-|,t-|)}). (3) 

We will show in the next section that the statistical function is related to the imaginary part of the Wigner-transformed 
Feynman propagator as follows : 

Gjic(w,k, x, t) — iImG_F(o;,k, x, t). (4) 

Thus the statistical function is purely imaginary. In literature, this is often defined without a —i factor in front as in 
. In that case it is real, and when integrated over uj at fixed k, x and t, it gives the time-dependent phase-space 
distribution of quasiparticles in the semiclassical limit at weak coupling. Here we will keep the —i factor - in this case 
it is also known as the Keldysh propagator. 

The statistical function/ Keldysh propagator can be indirectly measured by following space-time evolution of 
expectation values of operators, like conserved currents and the energy-momentum tensor. We note that both the 
non-equilibrium statistical and spectral functions need regularization for removing ultraviolet divergences, which are 
however expected to be independent of the state. What is relevant for physical observations is their dependence on 
the non-equilibrium hydrodynamic and non-hydrodynamic variables characterizing the state. 

At equilibrium, the statistical function is again independent of x and t due to translational invariance. Also it is 
determined from the spectral function using analyticity. The relation between the spectral function and the statistical 
function at equilibrium is known as the fluctuation-dissipation relation, which at thermal equilibrium is as follows in 
the bosonic case: 

G K (uj,k) = -i(n S E(v) + -).A(w,k). (5) 

It is called the fluctuation-dissipation relation because it implies the fluctuation-dissipation theorems like that 
relating the thermal electrical noise (current fluctuations) to the electric resistance. 

Away from equilibrium, analyticity is insufficient to relate the spectral and statistical functions. They evolve in a 
coupled way given by the so-called Kadanoff-Baym equations which can be derived from the two-particle irreducible 
action (see section II). In the semi-classical limit at weak coupling, the Kadanoff-Baym equations reduce to Boltzmann 
equation with quantum corrections. Formally these equations relating off-shell quantities are valid at strong coupling 
also. Even at weak coupling however, non-perturbative methods are needed to make complete sense of these equations 
and for making systematic expansions as mentioned earlier. In general we do not expect any generalization of the 
flucutation-dissipation relation away from equilibrium which can be stated in a universal way independent of the 
details of the non-equilibrium state and the nature of external perturbations. 

Numerical simulations of Kadanoff-Baym equations based on two-particle irreducible (2PI) effective action at weak 
coupling indeed show apparent thermalization - at long time the spectral and statistical functions satisfy Such 
simulations, as for instance, in O(N) scalar field theories at large N and weak coupling pij] , indeed show thermal- 
ization; but far away from equilibrium, there seems to be no simple fluctuation-dissipation relation between spectral 
and statistical functions. 

Summary of results : In this paper we study the non-equilibrium spectral and statistical functions in states 
which are perturbatively connected to thermal equilibrium with external sources absent. These are non-equilibrium 
states undergoing strongly coupled hydrodynamic and non-hydrodynamic relaxation. 

At strong coupling we expect a few operators to suffice in describing typical states. This is so because most of the 
operators except the relevant order parameters, conserved currents and the energy-momentum tensor will be expected 
to have large anomalous dimensions. This will be more generic in non-supersymmetric theories |L2], |l3| where we do 
not have chiral primary operators 3 , unless there are emergent symmetries at large N (as in [H]). 



The chiral pimary operators are special because they are typically protected from receiving large quantum corrections to their anomalous 
dimensions 
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Thus we can describe the non-equilibrium states holographically at strong coupling and large N using solutions 
using Einstein's gravity coupled to a few gauge fields (dual to conserved currents) and scalar fields (dual to relevant 
order parameters). Here we will further specialize to states where conserved currents and order parameters vanish in 
the course of evolution. Thus we will be able to describe these states using the hydrodynamic and non-hydrodynamic 
quasinormal modes of pure gravity along with their non-linear evolution. 

The holographic prescription for obtaining the spectral function in these class of states systematically in the deriva- 
tive (hydrodynamic) and amplitude (non-hydrodynamic) expansions has been obtained before. In this paper we will 
generalize this methodology to obtain the statistical function in the same expansions systematically. We will find that 
the non-equilibrium statistical function, like the non-equilibrium spectral function individually carries detailed imprint 
of the non-equilibrium state, particularly detailed information about the dispersion relation of the relaxation modes 
and their non-linear interactions. However, there is a simple generalized non-equilibrium fluctuation-dissipation rela- 
tion between the spectral and statistical functions, again set exactly by the final temperature of thermal equilibrium. 

This non-equilibrium fluctuation-dissipation relation being independent of the non-equilibrium state concerned is 
universal and is as follows : 

Gx(w,k,x,t) = -i(n B E(u) + -).4.(a;,k,x,t), (6) 

where tibe(w) is set by the final equilibrium temperature. We note the definition of both, the spectral and statistical 
functions involve Wigner transform, so this relation is non-local both in space and time. However, we recover the 
equilibrium fluctuation-dissipation relation at long times when the spectral and statistical functions are independent 
of x and t. Also this relation is strictly true in absence of external sources, so conservation of energy implies that the 
final equilibrium temperature is set by the total energy of the initial state. We elaborate later in section V that the 
relation is not in contradiction with locality and causality of the underlying field theory. 

Furthermore, at strong coupling we expect the Kadanoff-Baym equations to have strong quantum memory, thus 
a Boltzmann like limit where the spectral and statistical functions will behave quasi-locally is likely to be absent. 
Therefore, the non-equilibrium fluctuation-dissipation relation need not be specified by local data. What is surprising 
is that the only data required is the total conserved energy of the state in the large iV and strong coupling limit. 

This relation has been obtained in the probe-brane limit in , jl6| . 

Methodology in brief and organization of this paper : We describe briefly the methodology we will follow 
here along with the organization of this paper. 

We start from non-equilibrium field theory in section II. We review briefly the Schwinger-Keldysh formalism. We 
then obtain new parametrizations of the non-equilibrium fluctuation-dissipation relation which will be necessary 
for having consistent derivative/amplitude expansions of the spectral and statistical functions. We show that the 
parameters must satisfy additional field-theoretic constraints. Based on this non-equilibrium fluctuation dissipation 
relation, we find that any non-equilibrium Green's function can be written as an appropriate weighted sum of the 
non-equilibrium retarded and advanced Green's function. The weights are also subject to field-theoretic constraints. 
Additionally we obtain new formal understanding of some important features of the effective action for non-equilibrium 
Green's functions. 

In section III, we review the construction of gravity duals of non-equilibrium states and introduce the deriva- 
tive/amplitude expansions. We then study the dynamics of the scalar field dual to the bosonic order parameter in 
the same expansions. We find the unique boundary conditions for the non-equilibrium modes consistent with the 
derivative/amplitude expansions which give regularity at the horizon. The only new element here will be study of 
conjugation properties of the boundary conditions necessary for later analysis of field-theoretic constraints. 

In section IV, we map boundary conditions of the non-equilibrium modes to the parametrization of non-equilibrium 
Green's functions as an appropriate sum of the non-equilibrium retarded and advanced Green's functions, obtained 
in section II. 

This map is possible because of our previous work in which the boundary conditions for the non-equilibrium modes 
which give solutions regular at the horizon have been identified, and have been utilized to obtain the causal response 
function which gives the holographic non-equilibrium retarded Green's function |l3[. This builds on the holographic 
prescription for thermal retarded Green's function proposed by Son and Starinets [[L7j. We show the advanced response 
function obtained from the regular non-equilibrium solution gives the holographic non-equilibrium advanced Green's 
function which passes field-theoretic consistency tests. 

We then show that the sum of causal and advanced gravitational response functions evaluated using appropriate 
boundary conditions map to the parametrization of the non-equilibrium Feynman propagator required for it to have 
consistent derivative/amplitude expansions. The boundary conditions which depart from regularity are responsible 
for the non-equilibrium shifts in the weights of the non-equilibrium retarded and advanced Green's function in the 
weighted sum. We then study reality constraints on the parameters obtained from field-theory and show that the 
only boundary conditions which are consistent are those which give regularity at the horizon. Thus we find that there 
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are no non-equilibrium shifts in the weights of the non-equilibrium retarded and advanced Green's function in the 
weighted sum which gives the non-equilibrium Feynman propagator. 

We also give intuitive aguments why regularity at the horizon in combination with linear response theory should 
determine all non-equilibrium propagators holographically. 

In section V, we obtain the non-equilibrium statistical function from the non-equilibrium Feynman propagator, 
and hence the non-equilibrium fluctuation-dissipation relation. We argue why it's form is consistent with the locality 
and causality of the underlying field theory. We also give a conjecture on higher point non-equilibrium correlation 
functions in absence of external sources. 

In section VI, we show though our holographic prescriptions generalize to some extent in higher derivative gravity, 
we cannot expect our results to hold beyond strong coupling even when the classical gravity approximation is valid. 
The main reason for this will be that the bulk scalar will not be minimally coupled to gravity generically. 

In section VII, we conclude with a summary of the key results and a discussion on how the holographic non- 
equilibrium fluctuation-dissipation result can be tested in heavy ion collisions particularly. 

The appendices give relevant supporting material. 

Before proceeding further, we would like to mention that an interesting formalism has been proposed by van Rees 
and Skenderis in llq] for doing real-time holography by construction of a gravitational analogue of the Schwinger- 
Keldysh closed time contour. Though this formalism seems ideally suited for doing thermal field theory in real time, 
it is very difficult to adapt this to non-equilibrium states. The main difficulty lies in construction of smooth Euclidean 
caps for non-equilibrium geometries, which is necessary in this formalism to set non-trivial initial boundary conditions, 
i.e. to represent creation of non-trivial states holographically. 



Furthermore, there is an interesting approach [h9[ due to Son and Herzog for obtaining the thermal Schwinger- 
Keldysh propagators. It is difficult to generalize this to non-equilibrium context as well as this will involve construction 
of Penrose diagrams of non-equilibrium geometries. Indeed it is not known how to do this for non-equilibrium 
geometries in perturbation theory. 

However, elements of these approaches may be concretely realized in special non-equilibrium geometries. It will be 
also necessary to understand how to construct the non-equilibrium spectral and statistical functions holographically 
in presence of external sources using systematic expansions. We leave such investigations to future work. 

Notations : We will denote the position and momentum four-vectors without bold fonts, as for instance by x and 
k respectively. The spatial position and momentum vectors will be denoted in bold fonts, as for instance by x and k 
respectively. 



In this section, we will first briefly review the Schwinger-Keldysh formalism of non-equilibrium field theory. In 
particular we will define the effective action which gives the equation of motion of the expectation value of operators 
and their Green's functions. Then we will show formally this effective action contains no new information from the 
usual effective action, which is the Legendre transform of the generating functional of connected vacuum correla- 
tion functions. Finally we will establish some new relations between non-equilibrium Green's functions of bosonic 
Hermitean local operators. 



The basic tool of non-equilibrium field theory is construction of an effective action which is a functional of the 
operator and correlation functions. The evolution of the expectation value of the operator and the Green's functions 
in all states, including those which are away from equilibrium, can be obtained by extremizing the effective action 
T[0(x), G(x, y)}. By constuction, this effective action does not depend on any equilibrium or non-equilibrium variables 
like temperature, velocity or shear-stress tensor. These non-equilibrium variables parametrize the solutions, which 
can be obtained from definite initial conditions. 4 

The effective action T[0(x),G(x,y)] can be defined as the double Legendre transform of a generalization of the 
generating functional of vacuum correlation functions - Z[J(x), K(x, y)]. This generalized partition function includes 



4 Typically we need infinite variables to parametrize non-equilibrium states, hence non-equilibrium expectation values of operators and 
their correlation functions. However, we will see in the next section, we may expect that at least in field theories admitting holographic 
gravity duals, there are non-equilibrium states which can be parametrized by hydrodynamic variables and shear-stress tensor alone. 




II. ASPECTS OF NON-EQUILIBRIUM QUANTUM FIELD THEORY 



A. The non-equilibrium propagators from an effective action 
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not only a source J(x), but also an arbitrary bi-local interaction term K(x,y), and is defined as below : 
Z[J,K] = e lW[J > K] 



2?$ s exp 



+ J d 4 xJ(x)0(x) + ^ J d 4 xd 4 yO(x)K(x,y)0(y) 



(7) 



Above $ collectively denotes all the elementary fields. The operator O which couples to J locally and K bi-locally is 
a polynomial of these elementary fields $ and their derivatives. 

We then define the expectation value of the operator 0(x) and it's Green's function G(x,y) through : 

8W[J,K] 



5J(x) 
8W\J, K] 



6K(x,y) 2 



0(x)0(y)+G(x,y)\. (8) 



Eliminating J(x) and K(x,y) in favor of 0(x) and G(x, y), we can now do double Legendre transform to define the 
effective action as below : 

T[0(x),G(x,y)} = W[J,K}- j d 4 x J(x)0(x) - \ f d 4 xd 4 y K(x,y)(o(x)0{y) + G(x,y)\ . (9) 



J{x)- / d 4 y K(x,y)0(y), 



Clearly, 

STp, G] 
50(x) 

6T[0,G] 1 . , 

Therefore, in absence of sources J(x) and K(x,y), extremizing the generalized effective action T[0(x),G(x,y)] gives 
the dynamics of both the expectation value of the operator and their Green's functions. 

There is one important point in the above construction which we mention now. As stated in the introduction, for 
non-equilibrium states we need to know separately the spectral function (the commutator/anti-commutator) and the 
statistical function (the anti-commutator/commutator) for both bosonic/fermionic operators. 

To get complete information about the non-equilibrium Green's functions, we need to construct the partition 
function W[J(x), K(x, y)] and consequently the effective action T[0(x), G(x, y)\ over the so-called Schwinger-Keldysh 
closed real time contour C s , as shown in fig. [I]. This time contour travels from — oo to oo and then back from oo to 
— oo, thus transversing the entire real line first forward and then backward. 

In fact, the full closed-time contour ordered Green's function Gc 3 (x,y) can be written as a combination of the 
commutator and the anti-commutator, and hence as a sum over the spectral function and the statistical function as 
below 

G Ct {x,v) = GK(x,y) - ^A(x,y), sign Ca (x° - y°), (II) 

where sign Cs (x° — y ) is defined according to the contour-ordering over the full closed-time contour C s 5 . 

The usefulness of the effective action T[0(x), G(x, y)] is realized only when we can reformulate it diagramatically. 
This has been done first in H via a two-particle irreducible (2PI) effective action when the operator O is an elementary 
field. A full discussion of this is outside the scope of the present paper. We do mention however, that though there 
is a diagrammatic formulation, this is non-perturbative as the lines of the diagrams are the full propagator G(x,y). 
Therefore, the effective action T[0{x), G(x, y)] typically admits no perturbative expansion in a small parameter. 

Despite this drawback, in happy circumstances, we do get good approximations to experimental results or simula- 



When both the operator insertions are on the upper contour we get the Feynman propagator, for instance. In general both the 
commutator and anti-commutator are needed to get the propagator over the full contour. Thus the contour is convenient for obtaining 
the independent dynamics of the non-equilibrium spectral and statistical functions. 
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C s 



0(x±) 



DO 



— OO 



t 



0(x 2 ) 



FIG. 1: The closed time Schwinger-Keldysh contour C a is as above. The forward and backward directed parts of the contour 
have been displaced slightly above and below the real axis respectively just to distinguish them clearly. The operator insertions 
can be either in the forward or backward leg of the contour. 

tions, particularly when the coupling is weak and we can have a good intuition about the relevant diagrams to sum 
over in the effective action ||. We need to resum over infinte diagrams to remove infrared divergences. There is no 
systematic understanding about how we can improve a given approximation while removing the infrared divergences 
in order to better the description. Nevertheless, this is so far the best field-theoretic tool at hand which gives a 
fundamental description of a physical process in an approximation that is uniform in time, i.e. which not only works 
at initial time but also far in the future. Furthermore this also reproduces quantum corrections to the Boltzmann 
equation. 



We will like to demonstrate here that the effective action T[0(x),G(x,y)] formally does not contain any more 
information than the usual effective action T[0(x)}. We recall T[0(x),G(x,y)] is obtained from Z[J(x),K(x,y)] by 
definition as in eqs. (|sj) and @, while r[0(x)] is obtained from the usual partition function Z[J(x)], the generating 
functional of vacuum correlation functions. r[0(x)] is the usual one-particle irreducible (1PI) effective action when 
0{x) is an elementary field. The effective action r[0(x), G(x, y)] does not admit any perturbative expansion in a 
small parameter, but r[0(x)] does. 

The effective action T[0(x),G(x,y)] is necessary to obtain a description which is a uniform approximation in 
time, but our current understanding of such approximations seem to be inadequate. It is not known how to have 
an approximation scheme which simultaneously takes care of unitarity in a controlled manner, is renormalizable by 
standard counterterms and satisfies all Ward identities (for recent progress see pj). Therefore it is useful to have an 
argument for why in principle all the information is there in the usual effective action r[C(x)] for which such issues 
have been resolved. This will be useful for us later. 

To show this, it is sufficient to prove that W[J(x), K(x, y)} can be constructed from W[J(x)], as T[0(x), G(x, y)} 
and r[C(x)] are obtained from W[J(x), K(x,y)] and W[J(x)] respectively via Legendre transforms. From (||) we 
readily observe that 



B. Redundancy of information in the effective action 



SW[J,K] _ 1 / SW[J,K]5W[J,K' 
SK(x,y) ~2\ SJ(x) SJ(y) 



5 2 W[J,K] 
SJ(x)SJ(y) 



) 



(12) 



where we have used 0(x) — 5W[J, K)/5J(x) and G(x,y) — SO(x)/5J(y). Thus all functional K— derivatives of 
W[J, K] can be converted to functional J— derivatives of W[J, K]. Therefore W[J, K] can be constructed readily 
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from W[J, K = 0] = W[J]. This is also expected because indeed a quantum field theory is uniquely defined once we 
know all connected vacuum correlation functions, hence W[J]. 

C. New relations between non-equilibrium Green's functions 

We will now investigate the general structure of non-equilibrium Green's functions of local bosonic operators in 
non-equilibrium states. Here the operators will be taken in the Heisenberg representation. 

In particular, we will establish some useful relations between the retarded propagator and the Feynman propagator, 
based on the assumption that the non-equilibrium parts of both admit systematic perturbative expansions. 

The definition of the Feynman propagator is : 

g f {x u x 2 ) = -i(r (o{ Xl )o(x 2 )) y (13) 

where T denotes time ordering. 



We begin with the known relation (for proof see appendix A) 



(i 



1 

G F (x 1 ,x 2 ) = G K (x 1 ,x 2 ) - -A(xi,x 2 ) signal -t 2 ). (14) 



Note this is just a special case of (|ll|) when the operator insertion are in the forward part of the closed-time contour. 
However, 

- iA(xx,x 2 ) sign(t x — *2) = —iA(xi,x 2 )Q{ti—t 2 ) + iA(xi,x 2 )Q(t 2 —ti) 

= G R {xi,x 2 ) + Ga(xi,x 2 ). (15) 

Therefore, 

G F (x!,x 2 ) = Gjc(x 1 ,X2) + ]^{g r {x u x 2 ) + G A (xt, x 2 )j, (16) 
or in the Wigner transformed form 

G F {k,x) = G K {k,x) + ReG R {k,x), (17) 

using G R * (fc, x) = GA{k, x) which has been proved in appendix A. 

We further note that G/c{k,x) is purely imaginary which follows from it's definition as proved in appendix A. 
Therefore, we conclude from ( p"7j ) that 

ReG F {k,x) =ReG R (k,x), lmG F {k,x) = -iG K (k,x). (18) 

It is to be noted that unlike the retarded Green's function G R (k,x) or the advanced Green's function GA(k,x), 
the Feynman propagator is not analytic in w either in the upper-half plane or the lower-half plane. This can be 
seen from the instance of the vacuum Feynman propagator itself in which positive frequency modes are propagated 
in the future and negative frequency modes are propagated in the past, i.e. the vacuum Feynman propagator is 
Q(oj)G R (k) — <d(—uj)GA(k). The term Q(u)G R (k) implies lack of analyticity in cj in the lower-half plane as the 
positive frequency poles are shifted there with the standard it prescription, and similarly 0(— oj)GA(k) implies there 
is no analyticity in the upper-half plane as the negative frequency poles are shifted there. 

Generally speaking, if a function is analytic in a variable in the upper half plane or lower half plane, with other 
variables held fixed, it can be fully constructed from it's imaginary part or real part. For instance as ImG R (k, x) = 
— (l/2)A(k,x) (for proof see Appendix A), and as it is analytic in the upper-half plane in ui, 7 the retarded Green's 



6 From now on we will denote the commutator as A(x%, £2), where it is implicitly implied that the inverse Wigner transform has been 
done on the spectral function A(k,x). Similarly we will denote —i/2 times the anti-commutator as G/c(xi,X2). 

7 Analyticity of Gn(k, x) in lower-half u> plane is guaranteed by the 0(ti — (2) term in Gr(x\ , £2). 
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function can be constructed from the spectral function using 

G R {w,k,x,t) = - — — . (19) 

J_ 00 27T cj - w' + is 

In case of the Feynman propagator, it's real part Re G R (k, x) and it's imaginary part G;c(k,x) hold independent 
information in non-equilibrium states which cannot be kinematically constructed from the other because of non- 
analyticity in lj in the complex upper-half plane or the lower-half plane. 

We also note that as the Feynman propagator is symmetric, i.e. Gf[x\,x%) = Gp(x 2 , Xi), it follows that after 
Wigner transform that Gp(k,x) = Gp(—k,x) as proved in appendix A. 

The thermal Feynman propagator is a linear combination of the retarded or causal response with weight ube(w) + 1 
and advanced or anti-causal response with weight — tibe (w) for a given frequency, with neE(w) = J" - being the 

Bose-Einstein distribution at temperature T. Thus, 

G F (k) = (uheM + l)Gfl(k) - n BE (u)G A (k). (20) 

As T -> 0, 

n BE (cj) + 1 -> 0(w), -n BB (w) -> Q(-w), (21) 

so we recover the vacuum Feynman propagator in this limit. 

It is not hard to see that (£0|) satisfies Gp(k) = Gp(—k), because 

Ube(-w) + 1 = -tibeM, or - »be(-w) = »be(w) + 1, (22) 

and also Ga(-w, — k) = Gr(uj, k). 

We can rewrite the thermal Feynman propagator (pM) as : 



G F {k) = ReG R (k) +i(2n BE (u) + l\lmG R {k) (23) 
So, comparing (^) with the above we get 

G K {k) = ^(2r^BE(^) + l)lmG R (fc) (24) 

However, ImG fl (fc) = — (l/2)-4(fc). So, 

G K (k) = -i(n B v(u) + l)A(k). (25) 

The above is the fluctuation-dissipation relation in thermal equilibrium. It is so-called because it implies the 
fluctuation-dissipation theorems , like those relating the thermal electrical noise (current fluctuations) to the electric 
resistance. 

We will now try to parameterize non-equilibrium fluctuation-dissipation relation. Using this, we will get a param- 
eterization of non-equilibrium Green's functions. 

Let us assume we are near equilibrium. Then the non-equilibrium state can be parametrized in terms of hydro- 
dynamic and relaxational modes. This will allow us to do a perturbative expansion. Let us for instance linearize 
the state in terms of hydrodynamic fluctuations 8u and ST. In this approximation, both the spectral and statistical 
functions will depend linearly on these variables. Thus the non-equilibrium fluctuation-dissipation relation will be 
linear in this approximation. 

The exact nature of the dependence of the spectral and statistical functions on the hydrodynamic and non- 
hydrodynamic variables will be determined later in the holographic classical gravity approximation. As of now, 
we will take this dependence to be implicit. 

It is also clear we can improve the linear approximation by including quadratic dependence of the spectral and 
statistical functions on the hydrodynamic and non-hydrodynamic variables which we will do systematically later. The 
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general structure will be as follows : 

G K (k,x) = G ( ° q) (k) + G ( £ cq \k,x), A(k,x) =A ( - c ^{k)+A ( - nc(i '>(k,x), 
G { ™ q \k,x) = G^\k,x) + Gf{k,x) + ... + G < g\k,x) + ... , 

A^(k, x) = (k, x) + A^ (k, x) + .... + A^(k, x) + ... . (26) 

Above the superscript (n) denotes the order of multi-linear dependence of the term on the non-equilibrium variables. 
The non-equilibrium fluctuation dissipation relation can be parametrized as follows : 



(k) = -i(n BE M + ^)^ (e<l) (fc), 
(k,x) = -*[ d%f^°\k,k 1 ,x)A^\k 1 )-i(n BE {u;) + ^jA w {k,x),.... (27) 



Above f^'°'(k,ki,x) is real because each term in the expansion of A(k,x) is real while each term in the G/c(k,x) 
is purely imaginary. Also the Bose- Einstein distribution hbe(w) appearing above is always determined by the final 
equilibrium temperature. Furthermore as we have shown in appendix A, because Gfc{xi,x 2 ) is symmetric in x\ and 
x 2 , after Wigner transform it should satisfy Gtc(k,x) = G;c(—k,x). Similarly as shown in appendix A, A(k, x) = 
—A(—k,x). Therefore, these imply that 

fQ->°\k,ki,x) = f£'°\k,kt,x) +f£' 0) (k,ki,x), such that 
fg~'°\k,ki,x) and f^'° (k, k\, x) are real, and 
0) (fc, k lt x) = jf 0) (- fc, h ,x) , /i 1>0) (k, k lt x) = (-k,-k lt x). (28) 

Similarly, at the next order 

G%\k,x) = -i f d 4 k 1 d 4 k 2 f^ '°\k,k 1 ,k 2 ,x)A {cq \k 1 )A {cq \k 2 ) 

d i k 1 d A k 2 d A x 1 f {2X0) {k, ki, k 2 ,x, xi)^ (1) (fci, a;i)^ (cq) (k 2 ) 
-i J d 4 k 1 d 4 x 1 f^ 1) (k,k 1 ,x,x 1 )A {1) (k 1 ,x 1 )-i(n BE (uj) + ^jA {2) (k,x), etc. (29) 

Again by definition /( 2 ' 0: °'(/c, k\, k 2 , x), f( 2 > 1,0 '(k,ki,k2,x,xi) and j"*- 2 ' 1 - 1 (fc, k\, x, x\) have to be real and the Bose- 
Einstein distribution n BB (uj) appearing above is always determined by the final equilibrium temperature. Similarly 
the symmetry of the statistical function and anti-symmmetry of the spectral function in x\ and x 2 prior to Wigner 
transform give further restrictions on f( 2 ' -°\k, ki,k 2 ,x), /t 2 ' 1 * ) (fc, k\, k 2 ,x, x{) and / < - 2,1 - ) (k, k±,x, x\) as in d28j). The 
last two terms of ( p9| ) are linear while the remaining terms are non-linear corrections to the equilibrium fluctuation- 
dissipation relation. Clearly non-linear corrections are expected because the Kadanoff-Baym equations which give 
equations of motion of the spectral and statistical functions are non-linear. We easily note as G^ and A^ are 
quadratic in non-equilibrium variables, while G$ and A^ are linear, /( 2 '°>°) depend quadratically while /P' 1 ' ) and 
fi 2 ' 1 ') depend linearly on non-equilibrium variables. 

We will find the above parametrization of the non-equilibrium fluctuation-dissipation relation useful to determine 
it holographically. It can be also useful for experimental determination. 

It then follows from eqs. ([l8]) and (27) that the general structure of the non-equilibrium Feynman propagator for 
small perturbations away from equilibrium should be 

G F (k,x) = G { p\k) + G^ cq \k,x), G { ™ q \k,x)=G^\k,x)+Gf{k,x) + .... + G^(k,x) + ..., 



Gp q ^ 



(k) = (nBEH + l)G^ q) (fc)-n B E(^)G^ q) (fc), 



G { p\k,x) = iG < £ ) {k,x)+ReG ( n\k,x), 



d i k 1 f il ' 0) (k,k u x)(G { ° q \k 1 )-G ( 2 q) (k 1 )) 

(n BE H + l)G%\k,x) - n BE {uj)G ( l\k,x). (30) 
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d 4 ' 



We have also used A(k, x) — — 21m Ga(k, x) and Gr * (k, x) — G^fc, x) above. Similarly from ( p9| ) we can derive that 

Gf{k,x) = - f d%d%/C 2 »°.°>(Mi,fc,s)^^ 

%d A k 2 d 4 xifV< 1 >°\k,k 1 ,k 2 ,x,x 1 )(G%\k 1 ,x 1 )-G%\k u x 1 j) (G { ° q \k 2 ) - G ( ° q \k 2 )) 

+ j d*fci<i 4 n/( J ' 1 H* ) *i,i,ii)(4 1) (*i.n)-Gi 1 1) fe>n)) 

+ (nBE(^) + l)Gg ) (fc,x)-n B EMG^ ) (fc, a ;), etc. (31) 

We note that with the explicit form of fluctuation-dissipation relation (p7|), we can write any non-equilibrium 
Green's function as a functional of the non-equilibrium retarded Green's function. As the Wigner transformed 
advanced Green's function can be obtained from the Wigner transformed retarded Green's function simply by complex 
conjugation, an alternative representation is a weighted sum of the non-equilibrium retarded and advanced Green's 
function, where the weights are specified by the non-equilibrium fluctuation-dissipation relation. 

Thus any non-equilibrium Green's function G(k,x) should have the following structure in the linearized approxi- 
mation : 

G(k,x) = G (oq) (fc)+G (ncq) (fc,x), G (noq) (fc,x) = G < - 1) (k,x)+G (2 \k,x) + ... + G in) (k,x) + ..., 

G«(fc,x) = J d^f^HkM^JG^Hkt) + f i 2'°\k,k 1 ,x)G^\k 1 )) 

MG g } (k, x) + /£> (k)G 2> (k, x) , (32) 

where /^J^ should be determined by the equilibrium fluctuation-dissipation relation and Jr'a should be determined 
by in the non-equilibrium fluctuation-dissipation relation (p7|). Similarly, 

G^(k,x) = J d A k 1 d 4 k 2 ^'°\k,k 1 ,k 2 ,x)G^\k 1 )G { ^\k 2 ) 

+f^°\k,k 1 ,k 2 ,x)G^\k 1 )G ( r\k 2 )+f^ 
+/i^ 0) (fc,fc 1 ,fc 2 ,x)G^ ) (fc 1 )G^ ) (fc 2 )) 

+ J d 4 k 1 d 4 k 2 d 4 x 1 ^' \k,k 1 ,k 2 ,x,x 1 )G^\k 1 ,x 1 )G^' l) (k 2 ) 

+f^'°\k,k 1 ,k 2 ,x,x 1 )G^\k 1 ,x 1 )G { 2 q \k 2 ) + f^ \k^ 
+f i A m (k,k 1 ,k 2 ,x,x 1 )G { j\k 1 ,x 1 )G i 2 q \k 2 )) 

+ J d^k x d^x x (/jj 2,1,0) (fc, k u x, x 1 )G ( ^ > {k 1 ,x 1 ) + /i 2ll ' 0) (fc, k u x, x 1 )G { j\k 1 ,x 1 ) S j 

+/£ q) (lo)G^ (k, x) + /i cq) (k)Gf (k,x)), etc. (33) 

Above /jj 2 j?°'°\ e t c - should be determined /( 2 >°>°), etc. in (p9|). We will find the above form of the non-equilibrium 
Green's functions useful for their holographic determination. 



III. NON-EQUILIBRIUM GEOMETRIES AND THE REAL SCALAR FIELD 

In this section, we will first review the geometries which capture basic non-equilibrium processes like hydrodynamics 
and relaxation in the dual theory. Then we will study how the solution of a free massive real scalar field can be 
perturbatively obtained in these geometries. This real scalar field will be dual to a condensate, as for instance 
the chiral condensate of QCD. Finally we will determine which solutions of the bulk scalar field are regular at the 
horizon. We will show that the non-equilibrium corrections to the equilibrium solutions are unique, and determined 
by boundary conditions at the horizon which can be stated in a background independent manner, i.e. independently 
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of the details of the dual non-equilibrium state, in the regime of validity of perturbation theory. 

Most of this section is a review of the results of our previous works ^l| . The only new element that we 

introduce is a careful study of complex conjugation properties of various terms, leading to the background metric and 
solution of the bulk scalar field being real. 



A. Holographic description of non-equilibrium states 



The solutions of vacuum Einstein's gravity with negative cosmological constant, which settle down to black holes 
with regular future horizons, holographically depict many fundamental non-equilibrium processes. Each such solution 
maps to a non-equilibrium state via gauge/gravity duality at strong coupling and large N. Each such non-equilibrium 
state can be characterized by the expectation value of the energy-momentum tensor alone, because this determines 
the dual solution in gravity which is a regular perturbation of anti-de Sitter black brane completely J22| . 

Conceptually, these states are similar to certain special solutions of the Boltzmann equation called conservative 
solutions, which are completely determined by the energy- momentum tensor p0| . This is only a conceptual similarity 
because the Boltzmann equation is valid only in the weak coupling regime. Nevertheless, it is seen that one can 
systematically construct phenomenological equations for the energy-momentum tensor and any solution of these 
phcnomenological equations can be lifted to a full solution of the Boltzmann equation. Thus such conservative 
solutions give a systematic approach to obtain phenomenology of irreversible processes in terms of energy-momentum 
tensor alone. At long times such solutions can be approximated by purely hydrodynamic solutions of Boltzmann 
equation, known as normal solutions [ p3[ in literature. These normal solutions are thus special conservative solutions 
and allow us to compute linear and non-linear hydrodynamic transport coefficients at weak coupling and dilute 
densities, when the Boltzmann approximation is valid 8 . 

In the case of gravity, one can also systematically construct phenomenological equations involving the energy- 
momentum tensor alone, which when satisfied implies that the corresponding solution in gravity should have a regular 
future horizon 20 , 2l], |2f| . Just like in the case of the Boltzmann equation, these phenomenological equations can be 



expanded perturbatively in the derivative (hydrodynamic) and amplitude (non-hydrodynamic) expansion parameters 
as reviewed in appendix B. These equations have the most general structure valid in any field theory when the 
quantum corrections to the Boltzmann equation are considered. Also the values of phenomenological parameters are 
different at strong coupling. 

At long time, all regular perturbations of anti-de Sitter black branes become purely hydrodynamic. In such cases, 
the metric can be constructed by the now well-known fluid/gravity correspondence p6|-p8|. 

This hydrodynamic limit is also a natural consequence of the phenomenological equations which admit purely 
hydrodynamic solutions pp| , pl| , 25 . The phenomenological equations however also describe other non-equilibrium 



processes like decoherence and relaxation. 

Another special class of solutions to these phenomenological equations describe homogeneous relaxation plj . In 
this case, it has been explicitly proved that the future horizon in corresponding solutions in gravity is indeed regular 
as outlined in appendix B. 

In this paper, for illustrative purposes we will use two basic examples. The first example is that of the five- 
dimensional anti-de Sitter black brane perturbed by a hydrodynamic shear-mode. The corresponding metric is 9 : 



72 



2t {l - f (^) )5 Ui (k {h) y^r*e-^dUr 

■— 

'4r 



•73r( - V" *<W ^(k(h))e ik ^' x e-^* h (^) dzW) + 0(e 2 ), (34) 



where, 



/(s) = l-s 4 , h(s) = -ln(l-s 4 ). (35) 



Indeed such an approach is also valid for non-abeliangauge theories, like quantum chromo-dynamics (QCD) at temperature higher then 
the confining scale Aqcjj and low baryon densities fe4] . 
One can readily construct this metric using methods of |22| . 
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The first line in (|34|) corresponds to the AdS§ black brane at temperature T = ro/(7rZ 2 ), with I being the asymptotic 
radius of curvature and I 2 /r Q being the radius of the horizon. The second and third lines describe the hydrodynamic 
perturbations. Here e is the derivative expansion parameter which is k^/T. 
The corresponding energy-momentum tensor is : 

k 2 

V = diag(e,p,p,p)-i7?(A: (h)i <Ju(k (h ^^ (36) 

where 

Above e andp are the thermal energy-density and pressure, and r\ is the shear- viscosity satisfying i]/s = 1/47T (applying 
standard thermodynamic relation to determine the entropy density s). Also k 2 is given by k 2 = (87rGjv) -1 , where 
Gn is the five-dimensional Newton's constant. 

The conservation of the energy-momentum tensor is required for (34) to be a solution of vacuum Einstein's equations. 
In fact, this follows from the vector constraint of Einstein's equation. The conservation of energy- momentum tensor 
yields the Navier-Stokes' equation. The latter implies that the velocity perturbation (5u(k( h )) is transverse, i.e. 

5u(k (h) ) • k (h) = 0. (38) 

The above procedure can be extended non-linearly to yield higher derivative |27], ^8) and non-linear corrections |2§f| 
to the Navier-Stokes equations. 

The second example will be that of homogeneous relaxation. The energy-momentum tensor for homogeneous 
relaxation takes the following form : 

V = diag (e, p, p, p) + 7r^ nh) (t) . (39) 

It is easy to see that the above energy-momentum tensor is conserved for any 7Ty (t). In case of a conformal field 
theory, we further require 

^; h) w^ = o. (40) 

The above is also implied by the scalar constraint of Einstein's equation. 

It is convenient to do a Fourier transform over time to define 7r|" h ^ (w( n h)). The corresponding metric up to linear 
order in ?r l ( " h) (w (nh) ) is 

I 2 dr 2 I 2 t „ /rr n \ .n . n . o\ 2l 2 



dS = 2 



j^y + M - / + d - 2 + d v 2 ) + ^(-& h) (-(nh)) h (^Uo*}) + o{5% (4i) 



with 8 being the parameter of non-hydrodynamic amplitude expansion parameter - 7rj° /p. Furthermore, h(s, W( n h)) 



follows the equation of motion : 

d 2 M S , W(nh) ) ( 2 +(l + 3 ^> 3 - 6 ^ 4 ) dh(s,u {nh) ) oj 2 nh) l 4 



ds 2 s f( s ) ds r 2 \f 2 (s) 



\h(s,LJ {nh) ) =0, (42) 



such that 



/i(s,w (nh )) = s 3 + Q{s 4 ) as s -> 0. (43) 



This uniquely defines h. 

It can be shown that the ( pd| ) is regular at the horizon provided 7r|^ (k>( n h)) nas simple poles in uj^^ exactly at the 
location of the homogeneous quasinormal modes (^. These quasinormal frequencies have both real and imaginary 
parts, and unlike hydrodynamic modes both of these are large, i.e. 0(T). In case of the AdSs black brane these are 
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known to be |29j : 

w (n)± = ttT ± 1.2139 - 0.7775 i ± 2n(l T i) , for large n. (44) 

This follows from the general phenomenological equations (see apendix B). 

The above analysis can be readily extended non-linearly in the amplitude parameter - tt^^ /p, but we need to sum 
over all time derivatives at each order in the amplitude expansion to see manifest regularity (see appendix B). This 
is because the time derivatives of these modes are large, i.e. 0(T). This is usual for non-hydrodynamic relaxation 
modes even close to equilibrium. 

The general phenomenological equations as discussed in appendix B have covariant form. In our examples we have 
written them in the laboratory frame, where the final equilibrium configuration is at rest. This is useful from point 
of view of comparing our results for Green's functions with experiments. 

The derivative and amplitude expansions in the laboratory frame is as follows. The derivative expansion will count 
powers of k(b)/T, where k(b) will be the momentum of the background quasinormal mode. It will be useful also 
to consider the quasinormal modes off-shell (and impose the on-shell condition later). In that case, the derivative 
expansion will also count powers of ljoa/T, with lj^a being real and the frequency of the hydrodynamic quasinormal 
modes without the dispersion relation imposed. For a generic quasinormal mode, the off-shell frequency will be 
denoted by W(b)< The amplitude expansion will count powers of 8u, ST/T^ and 7Ty '/p( eq \ with all powers of 
ui(Y,)/T summed up at each order for non-hydrodynamic modes. 



B. The scalar field in non-equilibrium geometries 

Here we will study the a free real scalar field in the background non-equilibrium geometries described above. This 
real scalar field will be dual to an appropriate order parameter like the chiral condensate in the field theory. The typical 
temperature of the non-equilibirum state will be above that of symmetry restoration. Thus the chiral condensate will 
have no expectation value. Furthermore the chemical potentials will be assumed to be small, so the flavor and baryon 
currents can be neglected. These conditions indeed are valid for the quark-gluon plasma produced by heavy ion 
collisions at Relativistic Heavy Ion Collider (RHIC) in Brookhaven and at A Large Ion Collider Experiment (ALICE) 
in CERN. 

Thus our non-equilibrium backgrounds can be modeled by the regular perturbations of AdS black branes in vacuum 
Einstein's gravity. As we are interested in the two-point correlations in these non-equilibrium geometries, we can ignore 
the backreaction of the scalar field, which will contribute only to higher-point correlations. 

It is convenient to Fourier transform the dependence on the boundary coordinates. We first need to do this for the 
background metric. 

Let us take the example of the hydrodynamic shear-mode perturbation discussed before for concreteness. As 
evident from (|34|), to achieve this, we need to Fourier transform the velocity perturbation Su. We have already done 
the Fourier transform in the boundary spatial coordinates, but we haven't done the Fourier transform of the time 
dependence yet. 

The Fourier transform of Ju(x, t) can be defined as [ fl3"| : 

W (h) + *4^T 

such that (Su* (k( h ) ) = 5u(— k^))- If we integrate over the wn^ over the real line and then closing in with the semicircle 
at infinity in the lower half plane as shown in fig. 0, we find that 



dw (h) <5u(k (h) , W(h) ) e ^)- x e-^<") t = «5u(k (h) )e lk <M- x e-— (46) 

c 

Similarly in the case of the homogeneous perturbation given by ( |4l| ) we can write 7r|" (a;^)) as 

oo 

where d( n )ij are constants satisfying ai n )ij5ij = and W( n )± satisfies (p4|). 
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Re w (h) 



FIG. 2: Contour for integration C over uji h ) which picks the contribution from the pole in the lower half plane 



This can be generalized to any quasi-normal mode, which can also be represented as a solution of the linearized 
phenomenological equations. The dispersion relation of any quasi-normal mode can be written as : 

w (b)±(k( b )) = ±Rew (b) (k( b) ) -HIm<j (b) (k (b) ), Imw (b) (k (b) ) < 0, (48) 
This gives rise to a complex pole in the Fourier transformed 5u(w( b ), k( b )), 6T(lj^, k( b )), 7r l -" h ' ) (^(b)i k(b)) as below: 

1 \ <5u±(k (b) ) 



2iri J oj (h) ± Re w (b) (k (b) ) + ilm cj (b) (k (b) ) ' 
1 \ «±(k (b) ) 



5u(k( b ),W( b )) 

ST k(b) ' W(b) \2m) uj (h) ±Rew (b) (k (b) ) + Tmcj (b) (k (b) ) ' 

(nh),, \ _ -f_L\ a u±( k (b)) 

H I (b). (b)J l 2 7rJo; (b) ±Re W(b) (k (b) )+iIma; (b) (k (b) ) : 



aij±(k( b ))5y = 0. 



(49) 



We will denote the frequency and momentum of the background perturbation as W( b ) and k( b ) respectively to distin- 
guish them clearly from those of the scalar field which will be denoted by ui and k. Note both W( b ) and u are real by 
definition. 

The metric is real, and so are <5u(x, t), <5T(x, t) and tt^" (x, t). So we need to consider a pair of modes which arc 
complex conjugates to each other. The complex conjugation involves a subtlety which has not been mentioned in our 
earlier work. It involves 

• reversing the sign of k( b ) and W( b ) , as for example we can see readily from ( |45| ) that 



<$u(-k( h ),-w (h) ) 



= , j_x <m-V)) = rj_\ *u*(k (h) ) = Ju; 



(k(h),W( h )) 



(50) 



; (h) - l TPT 



• the integration over CcVm which makes the mode on-shell should be done over the contour C* as shown in fig. 
U which goes over the real line and closes in the semicircle over the upper half plane instead of the lower half 
plane, as for example 



c* 



du w 5u(-k(h), -W(h)) e - l V)- x e ^(h)t = <5u(-k (h) )e^ k ( h )- x e-^ t = <5u * (k (h) ) e - 4k < h >- x e - — *, (51) 



We note that the contour C* is a reflection of the original contour C over the real line. It is necessary to integrate 
the complex conjugates of (Esh over C* so that the complex conjugate modes also decay in the future. We readily see 
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Re w (h) 



FIG. 3: Contour for integration C* over cj( h ) which picks the contribution from the pole in the upper half plane 



that we need to combine the two modes to produce a real perturbation, as for example 

dw (h) <5u( W(b) ,k (b) )e lk ( b )- x e-^) t + f doj (h) 6u(-oj {h) , -k (b) ) e -* k < b > 'V^' 

= ^u(k( b ))e' ik<b) ' x e _ '' RG "( b ' (k ( b » )t + Ju(-k( b) )e _lk<b) ' x e lRc " (b)(k(b))t ^e Im " (b ' (k(b):it , 

for a general velocity perturbation 6u(x, t). 

The free scalar field satisfies the usual equation of motion in the non-equilibrium geometry which is 

(□ + m 2 )0(r,x,i) = 0. 



(52) 



(53) 



The d'Alembertian □ however is that of the non-equilibrium geometry. 

Let us first consider the case when the perturbation of the geometry away from the AdS black brane is linearized 
in Su, ST and ■ Let these background perturbation be a sum of the specific mode given by the frequency W( b ) 



and momentum — k 



(b)- 



and momentum k( b ) as in (|48|) and it's complex conjugate at frequency 

The profile of the scalar field can be constructed as follows. Let us consider a specific mode of the scalar field at a 
frequency uj and momentum k. Interacting with the quasinormal modes in the background, new modes at frequencies 
w±u;( b ) and momenta k±k( b ) will be generated. So in the linearized non-equilibrium background the general solution 
of the scalar field will take the form : 



(x,i,r) = A(w,k)($( oq )(w,k,r)e- iM - k ' x) + I dw (b) $ (ncq) (o;,k,w (b) ,k (b) ,r)e^ 4(( " + "( b ) ) ^ (k+k ' b ) ) - x) 



+ / ^ (b) $(" cq )( W ,k,- W(b) ,-k (b) ,r) e -' i «"-^ b )) t -( k - k < b ))- x ). 
Jc* 



(54) 



We note the following points regarding the above equation : 

• The non-equilibrium corrections given by <fr( neq \uj, k, W( b ), k( b ), r) and $( neq )(w, k, — W( b ), — k( b ), r) can be eval- 
uated systematically in the derivative and amplitude expansions. 

• Both $( neq )(w, k, W( b ), k( b ), r) and $( neq )(u;, k, — W( b ), — k( b ), r) will be determined uniquely by the boundary 
conditions at the horizon. 

• We need to perform integrations over u)n,) using the contour prescriptions described above so that the background 
is on-shell. 

• The most general solution can be obtained as above by simply superposing over different equilibrium modes at 
w, k along with their unique non-equilibrium corrections, with an overall arbitrary coefficient A(uj, k). 
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As an explicit illustration let us examine the profile of the scalar field in the background metric (|34|) with a 
hydrodynamic shear-wave perturbation. 

The equation of motion of the equilibrium part <&( eqL '(u>, k) is 



( u t?,w> + m 2 )<5(w' - uj)6 2 (k' - k)$( cq )(a;,k, r) 



= 0, 



where ^ S is the Laplacian in the unperturbed AdS black brane metric given by : 



+ ?''(?) 



d r + r 2 



to 



(55) 



(56) 



Above, / is the blackening function of the AdS Reissner-Nordstorm black brane which vanishes at the horizon located 
at r = I 2 /tq, and is as given in (p5|). 

The equations of motion of the non-equilibrium parts of the solution up to first order in the derivative expansion 
are : 

+ m 2 ) 5(u' -co- u {h) )S 3 (k' - k - k (h) )<i>( noq > ( W) k, W(h) , k (h) , r) = k, W(h) , k (h) , r) $< cq > ( w , k, r), 

+ m 2 )$(a/ - w + uj {h) )5 3 (k' k + k (h) )$( noq V, k, - W(h) , -k (h) , r) = V(oj, k, - U(h)l -k (h) ,r) 

$ (cq) (w,k,r), 

with 



(57) 



V(w,k,W( h ),k( h ),r) = Vi(w,k,a; {h) ,k {h) ,r) + V 2 (w, k, w (h) , k (h) , r) + 0(e 2 ), 

2r 2 

Vi(w,k,a;( h) ,k (h) ,r) 



9 

y 2 (^,k,w (h ),k( h ),r) = i— fe^^Jfe%A?(h)i<yuj(w(h)»k(h))* 



(58) 



Above, h(rr^/l 2 ) gives the hydrodynamic correction to the background metric which is proportional to knAiSuj + (if> 
j), and is as given in (|35|). Also has been expanded up to first order in the derivative expansion only, the first 
term V\ depicting the contribution at the zeroth order in the derivative expansion, and V% depicting the contribution 
at the first order in the derivative expansion respectively. Both are linear in 8\x as we have treated the background 
in the quasi-normal mode approximation. Note the derivative expansion counts only the order of appearance of 
hydrodynamic frequency ton^ and kp^ through derivatives of <5u, while dependence in lu and k should be treated 
exactly at each order. 

As indicated in (Q) , we need to integate over W( b ) so that the background metric is on-shell. After this integration, 
the profile of the scalar field is : 



,i(k+k (h) )-x 



*(x,t,r) = $ (cq) (w,k,r)e- 4( "*- k ' x) + ($ (ncq) (w,k,k (h) ,r)e-^ t , 

k 2 

+$ (ncq) (w, k, -k (h) , r ) e -^* e 4 ( k - k ("))- x ) e"^K 
Above, $( ncq '(w,k,k (h) ,r) and $( ncq )(w, k, -k (h) , r) are the residues of & ne <d(u>, k, u (h) , k (h) , 



5>( neq )(u;, k, — W(h), — k( h ), r) at the poles u) 



(h) 



Tik 2 h) /(47rT) respectively. $( noq )(w,k, 



(h; 



(59) 



r) and 
r) and 



$( ncq )(oj, k, — W(h), — k(h), r) depend linearly on Su(w^, k(h)) and 5u(— W(h)» — k (h)) respectively through the V 
terms in (57). The contours C and C* pick up the residues of the poles in 5u(cd(M, kp^) and <5u(— wp^, — kp^) in the 
upper and lower half planes respectively. 

It follows that <I>( ncq ) (w, k, k(M, r) takes the following form up to first order in the derivative expansion : 



$( noq ) (w, k, k (h) , r) = $i noq) (w, k, r) <5u(k (h) ) • k + ^ ncq) ( w , k, r) hkj k {h)i 5uj (k (h) ) + 0(e 2 ). 



(ncq) , 



(60) 



The similar expansion for $( neq )(w, k, — kp^, r) is obtained simply by reversing the sign of kp^ above. 

The above procedure can be readily generalized when the background metric is expanded non-linearly in the 
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perturbations parametrized by <5u, 5T and 7T ? -° 

Let us focus on the case of the non-linear hydrodynamic background whose explicit forms can be found in the 
literature. Let us consider the Laplacian when we take into account quadratic dependence on two distinct velocity 
perturbations Su(k^) and Su(k'^) for instance, at a given order in the derivative expansion m (i.e. at the mth order 
in the hydrodynamic momentum). The solution for $ will receive a correction quadratic in the amplitude of velocity 
perturbation, and at mth order in the derivative expansion this takes the form 

$< 2 ' m > (r, k, fc (h) , fc' (h) ) e^W+tfw)- + (fc (h) -> -fc (h) ) + (fc[ h) -> -k{ h) ) + (fc (h) -> -k [h) , k{ h) -+ -k{ h) ). (61) 
The radial dependence above can be determined the equation of motion : 

a£ BB 5 3 (k' - k - fc (h) - k[ h) )& 2 < m \r, k', fc (h) , k{ h) ) = S^ m \r, k, fc (h) , fc' (h) ), (62) 

where 0^ BB is the Laplacian for a scalar with four-momentum k in the unperturbed AdS$ black brane and 5'( 2 » m ) is 
the source term determined by the background perturbation. For m = 1 i.e. at first order in the derivative expansion, 
the source term S^ 2,1 ) can contain terms like (k • du(fc(h)))(k(h) • <5u(fc^))<I>( eq ), (k( h ) • Su(k'^))^ 10 ^ etc. Finally we 
need to integrate over uj^ and uj'^ over C or C so that the background metric is taken on-shell. 

Thus perturbatively the solution of the scalar field in an arbitrary background can be systematically expanded both 
in the derivative and amplitude expansions. 



C. The boundary conditions for regularity 

The behavior of the general solution $( cq )(tj, k, r) in the equilibrium black brane background can be split into one 
that is incoming at the horizon and another that is outgoing at the horizon. The incoming mode $ m ( eq )(u;, k, r) can 

.2 

be defined uniquely by it's behavior near the horizon r — ^ as below : 

^)(u,,k,r) „ (l-^)" 4 * (63) 
The outgoing mode $ out ( eq )(w, k) can be similarly defined uniquely by it's behavior near the horizon given by : 

*-W( W> k,r) « (l-^y* (64) 

Also, 

$ in (° q > * (w,k,r) = $ in (° q )(-w,-k,r) = $ out ( cq )( w , k, r) = $ out ( cq ) * (- w ,-k,r). (65) 

In fact, the effect of reversing the sign of k is trivial as the solution can only depend on it's modulus due to the 
rotational symmetry of the black brane background. 

It is natural to include only the incoming solution when u) > 0, because it respects the causal structure of the 
black brane which forbids anything coming out of the horizon classically. We can also follow |3(j] to argue that if we 
keep both incoming and outgoing modes, it will cause a singular back-reaction at the horizon. As we will be consider 
quasinormal mode perturbations of the background metric which are incoming at the horizon, we should exclude the 
outgoing solution of the scalar field to prevent singular backreaction at the horizon. Thus the general solution at the 
zeroth order is 

$( cq )(w,k,r) = A in ( cq )(w,k)$ in ( cq )(w,k,r). (66) 

Until the end of this subsection we will assume uj > 0. Later we will add a complex conjugate term necessary to 
make the full solution real. This complex conjugate will take care of the case w < 0. 

In order to study the regularity of the non-equilibrium solution let us first examine the approximation where 
the background metric perturbation is linearized in Su, 5T and 7r,-" h ' ) • In particular let us study the case of the 
hydrodynamic shear wave perturbation. 
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It is easy to see from (|57j) that the general form of the non-equilibrium part of the solution can be written as 

$( ncq )(w,k, W(h) ,k (h) ,r) = A in ^ ( W ,k, W(h) ,k (h) )<f in ^)( W + W(h) ,k + k (h) ,r) 

+ A out(nc q) ^ k W(h) ^ k(h) ) $ i„(c q) {u + W(h) ^ k + k(h) ^ r) 

+A in ^)( W ,k)$ in ( nc ^( W ,k, W(h) ,k (h) ,r) + ( W(W -w (h) ,k (h) -> - W(h) ). (67) 



up to first order in the derivative expansion. The first two lines above indicate the homogeneous solutions of (|57 
but now with oj replaced by u> + u>n^ and k replaced by k + k( h ). The coefhcients of these homogeneous solu- 
tions A m ' ncq )(w,W( h ),k, k(h)) and A out ( ncq )(u;, W(h), k, k( h )) respectively have to be linear in 5u and have a consistent 
derivative expansion also. Their dependence on ujn^ and k( h ) can be expanded systematically in terms of rotationally 
invariant scalars like Su ■ k, kikjkn^iSuj, uinAkikjkaAiSuj , etc. Up to first order in the derivative expansion, only the 
first two scalars will apear. The coefficients of these scalars should be functions of to and k only, as the depenedence 
on conA and k(M can be absorbed in coefficients of the scalars appearing at higher orders in the derivative expansion. 
Thus, up to first order in derivative expansion, we should have : 

A in ( ncq V, W(h) ,k,k (h) ) = 4 n(neq V,k) ,5u( W(h) ,k (h) )-k + 4 n(ncq) ( w ,k) fc i %fc (h) ^( W(h) ,k (h) ) + 0( £ 2 ), 
A out(„e q)(W;W(h)jk;k(h)) = A °ut(„e q)( ^ k) k(w(li)ik(|i)) . k + i ™'W (Wjk) fc l fc J fc (h)l <5 Mj ( W(h) ,k (h) ) + 0(e 2 )(68) 

We recall for the hydrodynamic shear mode Su- k.n^ = 0, so there are no more possible terms up to first order in kp^. 

Above 4 n(noq) ( W ,k^4 n(noq) (^,k), A out(noq) (w,k) and A out(noq) (w, k) are arbitrary. 

The last line of (|67j) denote the particular solution determined completely by the perturbation of the background 
( p8| ) and the zeroth order solution ( p6| ) as appearing on the right hand side of (|57j). There is no new constant 
appearing here. $ ln ( ncq )(u;, k, £*>(h), k(h), r) is the particular solution sourced by $ ln ( eq ^(a;, k, r) and can be thought of 
as the non-equilibrium correction to the incoming mode. It's behavior at the horizon is given by : 

•( > f ttTPV ujSu(uj,k (h) ) ■ k / Tr-oV' 3 ^ 

$'"(- q ) W ,k, W(h) ,k (h) ,r « *2[ ro m I 1 7F • (69) 

\ r Q J (2a; + W(h))w( h ) y I 2 J 

We observe this particular solution is simply proportional to the equilibrium incoming mode at the leading order and 
hence is regular at the horizon. 

Just like in the case of equilibrium, we should discard the outgoing solution to prevent potentially harmful backre- 
action. Therefore, we should choose 

A°^)( w ,k,a; (h) ,k (h) ) =0, i.e. A° ut(ncq) (u, k) = 0, for i = 1,2, .. . (70) 

We recall that as in (|54j) we need to integrate over W(h) on the contour C given in fig. |2|. It is clear from the form 
of A ln ( ncq ^ (cj, k, W(h), k(h)) m ( p8| ) that the pole in the lower half plane of (5u(w( h ), k( h )) as given by ( fl5| ) will produce 
a divergence at the horizon whose leading term is : 



O0 



(71) 



This can be readily seen from the behavior of the incoming homogeneous solution <l> ln ( eq ) (u; + w/M,k + k(M, r) near 
the horizon as given by (|63|), substituting LJnA with it's on-shell value. To get rid of this divergence, we also need to 
put 

A in(neq) ^ k> W(h)) k(h) ) = , i.e. A? (neq) {to , k) = 0, for i = 1,2, .. . (72) 

Similarly we can write the general solution for <I>( ncq )(w, k, — W(h)j — k(h), r) in (|57j), as a sum of homogeneous solutions 
given by $ in ( oq )(u;-u; (h) , k-k (h) , r) and $ out ( cq )(w-w (h) , k-k (h) ,r) with coefficients A in ( ncq )(w, -w (h) ,k, -k (h) ) and 
A out ( ne<1 \uj , — W(h), k, — k(h)) respectively. These coefficients are arbitrary but have a consistent derivative expansion 
as in @. The particular solution is A in ( oq )(w, k)$ in ( ncq )(o;, k, -w (h) , -k (h) , r) with coefficient A in ( oq )(w, k) and 
<j)Out(ne q ) ^ k _ w ^ _ k ^ r j which is completely determined by the equilibrium solution and is regular at the horizon. 

Once again, we should put ^4° ut ( ncc i) to zero to prevent harmful backreaction. We should also A m ^ ncq ^ to zero as the 
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integration over contour C* (see fig. |3j) as in ( |54| ) will produce a leading divergence at the horizon of the form (|7l|). 

In a way these boundary condition for regularity is also a consequence of reality constraints (|78| ) which relates 
^jm(neq) £ Q ^out(ncq) n0 ^ e as th e background metric is real, for every background (u>(h) , k(b) ) quasinormal mode, 
we need to include the complex conjugate (— W(b)j — k(b)) quasinormal mode. Thus ( [78|) tells us if the coefficients of 
the outgoing homogeneous solutions are set to zero, so should be the coefficients of the ingoing homogeneous solutions. 

The particular solutions which remain are regular at the horizon. Therefore regularity at the horizon uniquely 
determines the non-equilibrium correction to the profile of the scalar field in the unperturbed black brane background. 
The full solution is : 



A(w,k)($( cq )(w,k,r)e- l ^*- k - x ) + J dw (h) $ in ( noq )(o;,k,w (h) ,k (h) ,r)e~ t « tJ+w ( h )) t ~( k+k ( h ))- x ) 

+ f ^ (h) $ in ^( W ,k,- W(h) ,~k (h) ,r)e^(^"^)) t -( k - k 0'))- x )). (73) 



The above boundary conditions apply to all quasinormal mode backgrounds. In each case, we need to put the 
coefficient of the outgoing mode at the horizon to zero to prevent harmful backreaction. Also the coefficient of the 
incoming mode should be put to zero, otherwise the W(b) integration will produce divergence at the horizon whose 
leading term will take the form : 



^+Re„ (b) (k (b)) ) l m „ (b)(k(b) ) 

4tvT / \ 4-kT 

rr \ / rro 

P \ I 2 



The above is divergent because for any quasinormal mode Ima;(b)(k(b)) < 0. 

The particular solution is regular at the horizons in all the instances checked so far. As for instance in case of the 
metric perturbed by the homogeneous relaxation mode in (pi]), the particular solution vanishes at the horizon due to 
factors like : 

-?)"(h.(i. 

This particular solution gives the unique non-equilibrium correction to the equilibrium solution. The full solution as 




in (73) has only one arbitrary constant which is A(w,k). 



The above boundary conditions are also valid when we take into account non-linear dependence of the background 
metric on Su, ST and 7r|" h ^ . 

Let us consider the case of the scalar field in the non-linear hydrodynamic background. In particular, let us consider 
quadratic dependence on Su at m-th order in the derivative expansion. This non-equilibrium correction is <I>( 2 ' m ) whose 
equation of motion is (^). Clearly the general solution of <J>( 2 ' m ) near the horizon can again be separated into two 
homogeneous pieces, the incoming and the outgoing modes, and a particular piece which has no arbitrary integration 
constant and is completely determined by the source term S^ 2,m \ In order to prevent harful backreaction, we should 
put the coefficient of the outgoing mode to zero. Also as discussed before, the integration over the hydrodynamic 
frequencies W(h) and uj'^ will produce a divergence at the horizon for the incoming mode, as for instance in the case 
above with dependence on two hydrodynamic shear wave background modes like : 



. i <•> ( h ) ( h ). 



(75) 



Obviously the coefficient of the incoming mode has to depend on £u(k(h)) and Ju(k'^) as required by the order in 
the perturbation expansion as in (|68|). The poles in uj^ and uj'^ as in ( |45| ) produce the above divergent behavior 
after integration over contour C or C*. In a consistent holographic set-up, the particular solution should be regular 
at the horizon. 

Thus the full solution of the scalar field in the non-equilibrium background is undertermined only up to an overall 
constant A(u>, k). This is beacuse the equilibrium solution $ m (cj, k) has a unique non-equilibrium correction which is 
regular at the horizon. 

We also note that for the full solution of <&(r, x, t) to be real we need to superimpose also $ m (r, — w, — k) = 
3> out (r, w,k) with $ m (r, w,k). Note the outgoing boundary condition is indeed natural for the negative frequency 
modes. As long as there is CPT invariance in the dual theory, the outgoing mode with negative frequency is the 
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natural CPT conjugate of the incoming mode with positive frequency. 

We also note that the non-equilibrium correction to $ m (r, — w, — k) is also complex conjugate to the non-equilibrium 
correction to <I> ln (r, to, k). In particular the correction to $ ln (r, — ui, — k) due to the background quasinormal modes 
are complex conjugates to the correction to <& ln (r, — to, — k) due to the complex conjugated quasinormal modes, i.e. 

<& in ( ne «) * (r.w.k, W(b) ,k (b) ,^ b) ,k' (b) , ...) = $-(^)(r, - w> -k, - W(b) , -k (h) ,-cj{ h) , -k' (b) , ...). (76) 

As in the background perturbations for every (w( b ),k( b )) mode there is the complex conjugate (— uin,), — k(b)) mode, 
the full non-equilibrium corrections of $ m (w,k) and $ m (— w,— k) are complex conjugates of each other. We can 
readily see that the complex conjugate is also regular at the horizon. 
Similarly 

$°ut(ncq) + (r , Wj k ^ j k(b) ) ^ ( ) _ _ } = $ ou t (ne q ) ^ _^ _ k) > ^ > > ( _ ^ (??) 

For the more general solution which is not necessarily regular at the horizon, we also need 

A in ^)( w ,k, W(b) ,k (b) ,^ b) ,k' (b) ,...) = A out ( nc ^(^,k,^ (b) ,-k (b) ,^ b) ,-k' (b) ,...). (78) 

for the solution to be real. We can readily see these imply specific relations between coefficients of the derivative 
expansion of both sides of the above equation. We can also derive similar requirements at higher orders in the 
amplitude expansion. 



IV. BOUNDARY CONDITIONS AND NON-EQUILIBRIUM RESPONSE FUNCTIONS 

In the previous section, we have studied the solution of a free scalar field in a non-equilibrium geometry. These 
geometries represents non-equilibrium states which are perturbatively connected to thermal equilibrium in the deriva- 
tive and amplitude expansions. Utilizing these expansions, we will now build a framework in which we can determine 
the non-equilibrium corrections to the Green's functions holographically. 



A. On the physical significance of the boundary conditions at the horizon 



We have seen in the previous section that the non-equilibrium modification of the solution of the scalar field is 
determined uniquely by regularity. We required the non-equilibrium modes generated by the background quasi- 
normal modes to be regular at the horizon. It amounted to putting both the boundary conditions at the horizon, 
namely setting the coefficient of the incoming and outgoing homogeneous solutions to zero. 

We may wonder what is the physical significance of putting both boundary conditions at the horizon in terms of 



the dual field theory. This has been partly explored in 13 



Firstly, we indeed expect that the external source will receive non-equilibrium medium modifications just like the 
expectation value of the operator. This is because the external source as seen by the quasiparticles should be screened 
by the non-equilibrium collective excitations of the medium. Thus the Dirichlet boundary condition is the wrong 
thing to use for the non-equilibrium bulk modes as this would imply no such screening happens EJ. 

Secondly, the quasi-particle dispersion relation is modified by the non-equilibrium collective excitations of the 
medium. We know for instance the mass of the quasiparticles receive thermal contributions. Taking this thermal 
mass into account is necessary for curing infrared divergences in thermal field theory. In non-equilibrium field theory 
the effective mass is space-time dependent, which is natural given that the temperature is space-time dependent as 
well. Furthermore, the effective mass can depend as well on the velocity and shear-stress perturbations also. There 
is no systematic way to obtain the dependence of the effective mass as a function of the background temperature, 
velocity and shear-stress perturbations in non-equilibrium field theory. 

Our holographic prescription gives a way to achieve this. The quasiparticle dispersion relation in equilibrum is 
obtained by solving for k/ cq )(k) at a given k, such that the external source satisfies J^(ui^(k), k) = 0, the 
latter being determined by the incoming boundary condition at the horizon. This can be readily generalized to non- 
equilibrium to obtain the space-time dependent shifts of the quasi-particle dispersion relations as parametrized by the 
background velocity, temperature and shear-stress perturbations Q. This is reviewed in appendix C. 

The shift in the effective mass should come from resummation of infra-red divergences. As the infra-red physics is 
given by the dynamics of the horizon holographically, this indeed justifies why both the boundary conditions for the 



22 



non-equilibrium modes should be applied at the horizon. Thus we determine the non-equilibrium modification of the 
source uniquely and hence the spacetime shifts of the quasiparticle dispersion relations. 

Building on this intuition, we can claim that the dependence of the non-equilibrium contributions of the Green's 
functions on the temperature, velocity and shear-stress perturbations characterizing the collective excitations should 
also come from resumming infra-red divergences. Hence, we can build on the hypothesis that the boundary conditions 
which will determine the non-equilibrium contributions to the propagators should be applied at the horizon. 

Thus we will examine general boundary conditions at the horizon which need not lead to a regular solution in 
the bulk. In the next subsection, we will study how the non-equilibrium corrections to the external source and the 
expectation value of the operator are parametrized by these boundary conditions. Then in the following subsection, 
we will use linear response theory and field theoretic consistency to determine the relevant boundary conditions for a 
given non-equilibrium Green's function. 

We will find eventually that regularity at the horizon determines all the non-equilibrium Green's functions in the 
holographic classical gravity approximation. In the final subsection we will return to non-equilibrium field theory to 
intuitively understand this. 

B. Non-equilibrium modifications of external source and expectation value 

For the moment, let us assume we have specified a pair of boundary conditions at the horizon which need not 
give a solution which is regular at the horizon. This means we have made specific choices of A ln ^ neq \uj, k, W(b), k(b)) 

and j4 ln ( neq )(w, k, W(b), k(b)), the coefficients of the two homogeneous solutions. Any such choice has to be consistent 
with derivative and amplitude expansions. As for instance in (|68|), this amounts to making arbitrary choices for 

4 n(noq) (w,k) and A° ut(noq) (w,k). 

After such choices are made, the entire solution is determined uniquely. We can readily determine the normalizable 
and non- normalizable modes from the solution by following it's behavior at the boundary. The full non- normalizable 
mode, which is sum of the equilibrium and non-equilibrium parts, is identified with the external source coupling to 
the dual operator. Let us denote this as J{x). The full normalizable mode is identified with the expectation value of 
the dual operator in the dual state, namely 0(x). 

Thus the equilibrium part of the source and expectation value of the operator can be obtained from the behavior 
of $( eq )(fc,r) near r = 0, as below : 

$(° q )(fc,r) a J^(k)r A - A +O icq) {k)r A . (79) 

Above A is the anomalous dimension of the dual operator and is related to the mass of the bulk scalar field by 

A = 2 + v / 4 + m 2 P. (80) 

The non-equilibrium corrections to the source and the expectation value of the operator in the dual non-equilibrium 
state can be obtained from the behavior of $( ncq )(fc, fc( b j, r) and $' noq )(fc, — fc(b), r) near r = 0, 

$( noq )(fc, ±fc (b) , r) « J^(k, ±k {h) y- A + 0< ncq )(fc, ±fc (b) )r A . (81) 

In case of the hydrodynamic shear mode background, for example, the non-equilibrium source and expectation value 
takes the form as can be seen from (|67|): 

(J, O) (ncq) (cj, k, ± W(h) , ±k (h) ) = A in (" cq ) ( W) k, ±c (h) , ±k (h) ) (J, 0) in(oq) ( W ± W(h) , k ± k (h) ) 

+A out(nc q) kj ±W(h) ( ±k(h)) (Jj o) in(«0 {uJ ± W(h) ; k ± k(h) ) 

( W) k) (J, 0) in(ncq) (U, k, ± W(h) , ±k (h) ). (82) 

Above {J,Oy n(cq) ' out(cq) (o; ± w (h) ,k ± k (h) ) are obtained from asymptotic expansions of $ in ( oq ). out ( cq )( w ± 
W(h),k± k(h),r). Similarly {J, 0) ln ^ ncq \uj, k, ±w(h), ±k(h)) are obtained from the asymptotic expansions of 
$ in ( ncq )(w,k,±w (h) ,±k (h) ,r). We note each of $ in ( cq )- out ( cq )( w ± W(h) , k ± k (h) , r) and $ in ( ncq )( w , k, ±« (h) , ±k (h) , r) 
are solutions of the radial equation, therefore have the same asymptotic expansions. 

Furthermore, using the specific form of the coefficients A ln ' neq ^ out ( nct i) as given in (p8|), it is easy to see that the 
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non-equilibrium modifications of the external source and the expectation value take the form : 

(J,0) (ncq) ( Wl i W(h) ,k,±k (ll) ) = (J,0) { ™ q) (u,k) 5u(±w {h) ,±k {h} ) -k 

±(J,O) { ™ q) (i0,k) k l k ] k {h)l Suj(±tu {h] ,±k {h) )+0(e 2 ), (83) 

with 

(J,ot cq \u,k) = 4 n(noq) (w,k)(^0) in(oq Vk) + Ar t(no %,k)(^,0) out(oq) ( w ,k) + (J,Op acq) (Lo 7 k). (84) 

In this form it is explicitly clear how the non-equilibrium modifications of the source and the expectation value depend 
on the boundary conditions at the horizon, namely in the choice of ^4™( ncq )' out( - ncq - ) (a;, k). 

The full space-time dependence of the source Jix) and the expectation value of the operator 0(x) thus take the 
form : 

(J,0)(u,,k,x,t) = (J,0)^( W ,k)e i < k -- wt > + [ duj {h) {J, k, W(b) , k(b) ) e i((k+k (b3 )-x-(^ (b) ) t ) 



+ / dw (b) (J 1 0)W( W) k,- ll , (b)l -k (b) )e i (( k - k ('))*-( w - w wW. (85) 
Jc* 

After the integrations over W( b ) required to take the background on-shell are done as above, we are left with the 
residues at the poles of the quasi-normal modes. Therefore, 

(J,(D)(u),k,x,t) = (J,0) {cq \uj,k)e l{k - x - ut) + (^ ;0 )(ncq)^ k ^ b) ) e4( (k+k (b) ^ 

+(J r ,0) (ncq) (w,k,-k (b) )e t(( ^ (86) 
Thus in the case of the hydrodynamic shear-mode background we get 

k 2 

(J,(D)(u,k,x,t) = ( l 7,G) (cq) (w,k)e ,; ( k - x - wt ) + ( l 7,0) (noq) (w,k,k (h) )e 4((k+k ( h ) ) - x)) e- l " t e-^* 

+(J,0) (ncq) ( Wl k,-k (h) ) e , « k - k w> x » e - ut e-^ t I (87) 

where 

(J,0) M ( U ,k,±k (h) ) = (J,0)[ ncq) (u,k) 5u(±k {h) )-k 

±{J,0) { ™ q) {u:,k) kikjk^Su^k^+Oie 2 ), (88) 

with {J, 0) • neq ^ given by (jS4|). Comparing (|S3| ) and ([§8|), we see that we have extracted the residue of 5u(±oj( h ), ±k( h )) 
as given in ( ff!i| ) at the poles W(h) = ±ik^ / (4ttT) . In general {J, 0) < " ncq \uj, k, ±k(b)) will take the form 

(J, O) (ncq) (w, k, ±k (b) ) = A hl ( ncq ) ( W , k, ±k (b) ) (J, 0t {cq) (w, k) + A° ut < ncq ) (w, k, ±k (b) )(J, 0) out(eq) (w, k) 

+ (^,0) in(ncq) ( W ,k,±k (b) ), (89) 

with A m ( ncq )' out ( n °q) and (J',0) ln ^ ncq> having consistent derivative and amplitude expansions. As for instance 
^ in (" cq )(^,k,k (b) ) =4 i ; (neq) (a; ) k)7r ii (k (b) )fc i ^ + ... . 

The above discussion is readily generalized by including non-linearities in the dynamics of the non-equilibrium per- 
turbations of the background. As for instance, if we take into account the quadratic dependence on two hydrodynamic 
shear modes <5u(k( h )) and 5u(k'^), the corrections will take the form, 

/ k (h) + k (h) A 

(^,0)("^( w ,k,k (h) ,k' (h) ) e l (V> + V))- x e — * e V ^ J + ( k(h) _>_ k(h) ) + (k' (h) ->-k' (h) ), (90) 
after doing integrations over and uj'^ . Furthermore, 

(J,0)^)(o;,k,k (ll) ,k' (h) ) = ( l 7,0)|" e9) (a;,k) (<Su(k (h) ) ■ k) (ju(k' (h) ) ■ k (h) ) + .... , (91) 
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at the first order in the derivative expansion. All these terms are fixed uniquely in terms of the boundary conditions 
at the horizon. 

The general form of the non-equilibrium corrections to source and expectation value as in ( |89| ) after including 
non-linear dynamics of the background involves 

{J, C) (ncq) ( w , k, ±k (b)1 , ±k (b) „) = A m ^{u, k, ±k (b)1 , ±k (b)n )G7, 0) m ^\uj, k) 

+A out(„e q)(w5 ^ ±]f(b)i) ^ ±k(b)n )( i7> 0)° Ut ^)( W , k) 

+ (J, Oy a ^(u, k, ±k (b)1 , ±k (b)n ), (92) 

which will appear with a spatio temporal factor 

e i(±k (b)1 ±...±k (b) „)-x e -i(±Rccj (6)1 (k (b)1 )±...±Rcw (b) „(k (b) „))t e (±Imw(i, )1 (k (b)1 )±...±ImLj ( i, ) „(k( b) „))t 



C. Mapping linear response to non-equilibrium Green's functions 

There is a basic linear response to an external perturbation in any physical system - the causal response. This 
causal response gives the retarded Green's function. 

Holographically, the equilibrium causal response maps to the incoming boundary condition at the horizon Jl7j . This 
is expected to be the case because it is the only boundary condition that is consistent with the causal structure of 
the black hole which forbids anything propagating out of the horizon classically. 

We have seen that the incoming equilibrium solution of the bulk scalar field has a unique non-equilibrium correction 
which is regular at the horizon and is free of potentially harmful back- reaction at the future horizon. We need regularity 
at the horizon to preserve the global causal structure of the black brane. We recall that this non-equilibrium solution 
which is regular at the horizon requires very specific boundary conditions at the horizon. These boundary conditions 
have been determined even after including non-linearities in the dynamics of the perturbations of the background 
quasi-normal modes. 

Furthermore, this non-equilibrium solution which is regular at the future horizon uniquely determines the non- 
equilibrium modifications to the external source and expectation value of the dual operator. It has been proposed 
that it determines the non-equilibrium corrections to the causal response, and therefore the non-equilibrium retarded 
Green's function | Q . 

The retarded Green's function in any state is given by the causal response in that state. Thus the holographic 
non-equilibrium retarded Green's function is : 

Gflte.iuxata) = / du^k ^"'** 1 '* 1 ] . (93) 
J J(w,k, X 2 ,*2) 

Let us go back to the linearized approximation when the non-equilibrium perturbation in the background is given 
by a single quasi-normal mode or a single relaxation mode in the dual theory. Then the full spatio-temporal form of 
J and O above is as in (^6|).The boundary conditions required for regularity at the horizon are that the coefficients of 
the homogeneous solutions j4 m ( ncc i)> °ut(ncq) should be set to zero, so that the non-equilibrium corrections to external 
source/ expectation value are given fully by (J, 0) m ^ ncq '(u;, k, ±k( b )) respectively. 

One can see that after doing the Wigner transform (for intermediate steps see appendix D) the non-equilibrium 
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retarded Green's function takes the following form [fU 



G fl (w,k,x,t) = j dun f <PkiG%°(wiM) 



<5(w-wi)(5 3 (k-ki) 



X(o; 1 ,k 1 ,k (b) )5 3 fk-k 1 - ^ . , 

kr 



-A^ 1 ,k 1 ,k (b) ),5 3 (k-k 1 + -^ 



ik (b )-x -iRew(b)(k {i , ) )t Imw (b) (k (b) )t 



+_L f M{uiM, -k(b))<5 3 fk - kj + ) ? -L 

27TI y W/ \ 2 / A , | Rc^(b)(k (t ,)) | . Im^ (b) (k (b) ) 



-^(Wi.kx.-k^^^k-k!-^ 



e -ik ( b) 'XgiRe u w (k (b) Jtglma^ (k (b) )t 



- Wl - Rc "(b)( k C')) _ ^ Im "(b)( k (b)) ^ 



(94) 



with 



R y ' J(°q)(w,k) v ; 

being the equilibrium retarded Green's function and 

^inCneq)/ k ±k \ jin(n Cq ) / k ±fc \ 

^k^H 0( ; q)( ; 1;kl) (b)j . ^ 1)kl ,k (b))= J( : q)( : likl) (b)j - (96) 

The complex conjugate mode corrections included above needed for complete field-theoretic consistency (see below) 
were not included in our earlier work [ pL~3| . 
The crucial elements to note are : 

• The non-equilibrium contributions involve two pieces, spatio-temporally comoving with the background quasi- 
normal modes at momenta ±k( b ) respectively. 

• The non-equilibrium contributions to the equilibrium retarded propagator at frquency uj and momentum k 
involves a convolution with support at momenta k ± k( b j/2 and poles at frequencies oj ± Reoj( b )(k( b ))/2 ± 
Tmw (b) (k (b) )/2. 

These features reflect the classical gravity approximation. We will see how these generalize beyond the quasi-normal 
mode approximation of the background. 

Also note that as the non-equilibrium correction to the equilibrium solution which is regular at the future horizon 
is unique, the overall solution is undetermined up to an overall multiplicative constant A(u>, k) as in ([54]). This mul- 
tiplicative constant A(ui,k) cancels between the numerator and denominator in (|93|). Therefore the non-equilibrium 
retarded Green's function is uniquely determined by this holographic prescription. 

The advanced Green's function should be naturally related to the advanced response. Thus, 

J J(w,k,xi,ii) 

The above also follows from the definitions that GU(xi, t\\ X2, tz) = Gr(x2, t%\ xi, ti) (for proof of the latter see 
appendix A). For the holographic definition of advanced response, both O and J above must be evaluated in the 
same non-equilibrium solution as in ([)3|) in the case of the retarded propagator, i.e. in the solution involving the 
unique correction to the equilibrium incoming mode which is regular at the future horizon. Therefore, O and J take 
the form as in fl8§) with (J,0) (ncq) = (J , 0) in(ncq) . 
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From appendix A, it follows that after the Wigner transform Ga(uj, k, x, t) — Gr(— to, — k, x, t). Therefore 
G^(a;,k, x, t) is given by with (w,k) replaced by (— u), — k). 

Also G,r * (fc, a;) = G^ {k, x) as proved in appendix A. This is also satisfied as can be seen from ( |9~4| ) using Ga (k,x) — 
Ga(—k,x). To do this one first reverses the sign of the integrated variables u)\ and ki and then uses results from 
section III.C which imply that (J,0) {cq) *(k) = (J,0) {cq) (~k) and (J, 0) in(ncq) *(fc, k (b) ) = (J, 0) in(ncq) (-fc, -k (b) ). 
Thus our holographic prescriptions for the non-equilibrium retarded and advanced propagators pass important field- 
theoretic conssitency tests. 

At the beginning of this section, we have already made the hypothesis that non-equilibrium Green's functions in the 
holographic classical gravity approximation should be determined by appropriate boundary conditions at the horizon. 
These boundary conditions should be determined by the consistency of the map with field-theoretic requirements. 

In order to map non-equilibrium Green's function to boundary conditions at the horizon, we start from the gravity 
side. Consider an arbitrary solution of the bulk scalar field in the non-equilibrium geometry which is consistent with 
the derivative/amplitude expansions but not necessarily regular at the horizon. This means that the space-time profile 
of the source and expectation value takes the form ( |86| ) with the non-equilibrium parts taking the general form ( |89| ) . 
We should make a choice of both A m ( neq )(u;, k, k( b )) and ^4 out ^ ncq ^(w, k, k( b )) to specify the solution uniquely. This 
amounts to making choices of coefficients of J 4^ n( - noq ' ) (aj, k) and ^4° ut ( ncq ) ( Wj k) Q f terms like 5u(k( b )) • k, TT^ h \]s.^)kikj 
in their derivative / amplitude expansions as in (^8|) . 

The most well-defined response functions on the gravity side are the causal and advanced response functions. Let 
us consider the following response function evaluated in an arbitrary non-equilibrium solution : 

ni \ ( j rfu( f{ e i)t xO(w,k,xi,ti) (oq) 0(w,k,x 2 ,t 2 )\ , nQ , 

G(x 1 ,x 2 ) = dujd*k\ f y R H '(uj)— — — + f A (uj)—, — r — H . (98) 

J V J(w,k,x 2 ,i 2 ) J(w,k,xi,ii) I 

Physically the boundary conditions given by choice of A m ( ncq )(u;, k, k( b )) and A ont ^ neq \uj, k, k( b )) can be interpreted 
as follows. The full response function above at equilibrium is the sum of the causal response with weig ht/jj eq) H and 
the advanced response with weight (w) . The boundary conditions ^4 m ( ncc i) (w, k, k( b j ) and ^4 out ( nct i) ( W) k, k/ b ) ) thus 
shift the weight of the non-equilibrium causal response and the non-equilibrium advanced response in the total weighted 
sum, because individually they are given by the boundary conditions A m ( ncq \uj, k, k( b ) ) and A out ( ncq >(w,k,k (b) ) set 
to zero. 

It is then easy to repeat the steps which lead from (|)3|) to (Q) as explicitly shown in appendix D. We find that ( |9^ ) 
reproduces our general parametrization of any non-equilibrium Green's function as a weighted sum of the retarded and 
advanced Green's functions (|32]), which is valid at the leading order in amplitude expansion, with the identification 

f(eq) 



Rcu (b) (k (t) ) „- Im"(b)(k (b )) 
W-Wi 5 1 2 



+A/> 1 ,k 1 ,k (b) )5 3 k-k 1 + ^ 



3/V , k (b)\ 1 



2 / ( , , , , Rc^ (b) (k (b) ) Im^ (b) (k (b) )) 
LO — U)\ -\ g r I 2 



g ik (b) -Xg-iRe u (b) (k (6) )t glm w (b) (k (b) )t 

•^M (m (u^M, -k (b) ) <5 3 (k - kl + ^) 



2;r/ \" ■ ' '"" l - 2 J ^ ^ , Rc^(k (t) ) | , im^ (b )(k (b) r 



+N ( Wl! k 1 ,-k (b) )5 3 (k-k 1 -^l 



2 ^ (co ul Rc "(b)^ k ('')) „- Im "(b)( k (b>) ^ y 

e -ik (b) XgiRe w (b) (k (i)) )i g Im w (b) (k (b) )t ^ (gg\ 
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where 



.M(wi,ki,±k (b) ) 



^in(ncq) : ki ^ ±k(b) ) -in(eq) ]f x ) + A out ( necl ) (_ X , kl , ±k (b) )£> out ( cc l) kl) 
= OM(^kO : 

k +k s ^ (ncq V, k 1; ±k (b) ) - J^Hun, ki, ±k (b) ) 
AA^.kx.ik^) - _____ 

^in(noq) ^ ^ ki ^ ±k(b) ) Jin(eg) ^ ^l)_h A out ^ (- 1 , kj , ±k (b) )^ out C_0 (_! , kj ) 

jCeq)( Wl| k!) 



(100) 



Also f^ {k,ki,x) is given by the above equations (|99| ) and ( |l00|) with fj^(u)) replaced by /^^(lu) and (w,k) 
replaced by (— cj, — k). 

Thus there is a map between between parameters /jj^J in the leading order parametrization ( [32] ) of non- 
equilibrium Green's functions required for consistent pertubation expansions and leading order boundary conditions 
^m(ncq), out(noq) £ or non . equilibrium modes set at the horizon. However, f^'^ should also satisfy other field- 
theoretic requirements, aside from producing consistent derivative and amplitude expansions. 

Let us work this out explicitly for the non-equilibrium Fcynman propagator. In this case we should set in the gravity 
response function ( |98| ) /j^ (w) — hbe^) + 1 and /ji (_) — f^ q \— _) = -)Jbe(w) to reproduce the equilibrium 
Feynman propagator. 

The holographic non-equilibrium Fcynman propagator Gf(x\, x 2 ) is therefore given by 



G F {xi,x 2 ) = J diod 3 k ((nB&(u) + 1 



0(_,k,xi,ti) 0(_,k,x 2 ,i 2 ) \ , ini , 

— — : — r-n BE (u))— — —r- , (101 

J (_,k,x 2 ,i 2 ) J(„,k,xi,ti) I 



with both J and O evaluated in a non-equilibrium solution with boundary conditions given by A 1D ^ neq \oj, k, k( b )) 
and A out ( noc i)(w,k,k (b) ), both of which we should determine now using field theoretic consistency conditions. 

Note the holographic Feynman propagator ( |101| ) is symmetric in x% and cc 2 as required. This can be seen by 
exchanging x\ and x 2 while reversing the signs of the integrated variables to and k. 

Let us now compare the holographic form of the non-equilibrium Feynman propagator ( |101[ ) with the parametriza- 
tion ( |30| ) that yields consistent perturbative expansions obtained earlier. Our comparison yields that 

f Mf kikux) = _< r _: P ( fn 0) ^ _ + __ ^ _ _ _ ) 

-5(/r(*.*i.*)-/]i ,0) (*.*i.4 ( 102 ) 

where fx\k, k\,x) is given by ( p9| ) with /jj (_) = iibe(w) + 1 and (k, k%x) is also given by ( p9| ) with /jj Cq ^(cj) 
replaced by /^ eq ^(u;) = /jj Cq ^(— w) = — ube(^)- Also P denotes the principal value. 

For the case of the Feynman propagator, using (frl ) and /Ji (w) = /jj 013 ^— — ) we find that 

fg>°\k,ki,x) = f£<°\-k,-ki,x). (103) 
We can now compare with the requirements ( |28| ) to obtain 

$'°\k,ki,x) = 0, 

f£'°\k,h,x) = fW{k,h,x). (104) 

We recall that the above implies that the constraint of symmetry in x\ and X2 is satisfied for the holographic 
Gf(xi,x 2 ), which is also evident in the defining form (101). 

There are two crucial field theoretic requirements in ( psj ) which should be additionally satisfied in the above 
identifcations. These arc that fg'°\k,ki,x) and fj^'°\k,ki,x) should be real. One can readily see from ( p9[ ) that it 
is impossible to satisfy these for all x unless both /j? 1,0 ^ (k, k%,x) and fjk' ^ (k, k\,x) both vanish identically. The latter 
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is possible only when A ln ( ncq,> and yl out ( nc q) are both set to zero giving the solution which is regular at the horizon 



This implies there are no non-equilibrium shifts in the weights of the causal and advanced responses in (101) 
Thus we have proved when the quasinormal approximation of the background is valid, 



G F (k,x)= (tibeM + l)G R (k,x) - n BE (uj)GA(k,x). (105) 



This simply follows from 101 using boundary conditions where j4 m ( nc< i) and ^4 out ( net i) are set to zero in ( |l~0l| ) and then 
doing Wigner transform. 

We can now go beyond the quasinormal mode approximation of the background. This means taking into account 
non-linearities in the dynamics of the hydrodynamic and relaxational modes constituting the non-equilibrium state. 



We note that (93) still gives the holographic retarded Green's function because it follows from linear response theory 



- the expectation value and source are obtained from the full non-equilibrium solution which depends non-lincarly on 
the non-equilibrium perturbations <5u, ST and Trf" 1 ^ and is regular at the future horizon. We have discussed how we 
obtain this full non-equilibrium solution in section III.C - namely we apply the same boundary conditions which cut 
down both the homogeneous incoming and outgoing solutions at each order in the amplitude perturbation expansion. 

The retarded response function depends non-linearly on each term in the perturbation expansion though. Let us 
denote (n) as the order in the amplitude expansion, i.e. it counts the total number of pertubation parameters <5u, ST 



and 7Ty as in section III.C. Then the retarded Green's function up to ?i-th order in the amplitude expansion is 
obtained from 

£)( c q) _|_ QW _|_ _ _|_ Q(n) 
j(eq) + J(l) + _ + J (n)- 

At the n-th order we not only get / ' but also {O^ / J^) n . Above we have suppressed the derivative 
expansion, but it clear we can do a double expansion in the derivatives and amplitudes of the background perturbation 
close to equilibrium. 



Similarly the advanced Green's function is given by (97), but the expectation value and the source are obtained 



from the full non-equilibrium solution which is regular at the horizon. One can readily check that the full holographic 
retarded and advanced propagators also satisfy the relation Gr * (k,x) = Gyi(fc,x) simply repeating the same steps 
discussed earlier to all orders in the perturbation expansion. 

One can now map the holographic linear response as given by (|98|) to a given non-equilibrium propagator. In d98| ) 
we need to use an appropriate boundary conditions at the horizon instead of those fixed by regularity at each order 
in the amplitude expansion. These boundary conditions give the non-equilibrium shifts of the weight of the retarded 
and advanced Green's functions which consitute the given non-equilibrium propagator. 

It is also easy to see that, at the second order in the amplitude expansion, we will obtain the same structure 
( |33| ) as we argued from field theory from the holographic response (|98|). Namely f^'^ A ar aa(^^ ^i, k 2 , x) will 
depend on the boundary conditions on the solutions <j>( 2 <™' m ) ; where n and m denote the order of the deriva- 
tive expansion acting on each of the two amplitude expansions. These are reflected by the boundary condi- 
tions given by A^ n ' m ^ n( - no ^(k,±k {h)1 ,±k {h)2 ) and A^ n ' m ^ out ( nc ^ {k, ±fc (b)1 , ±fc (b)2 ). We note that the reality con- 
straints (|78|) tells us there are exactly four independent boundary conditions for the four independent parameters 

fuR RA ar aa(^^ ^i' ^2 5 x). Without writing the detailed expression for /rr'ra ar aa(^^ ^2,%) we can readily see 
that taking all orders in derivative expansion into account, 

• the dependence on k is given by a sum of four terms which get supported by <5 3 (k — ki — k 2 ± k(b)i ± k(b) 2 ) 
respectively, 



as function of U), these terms have poles in u> — u>\ — u 2 ± w(/c(b)i) ± w (^(b)2) respectively, where w(fc(b)i) and 
w(fc(b)2) are the complex dispersion relations of the two quasinormal modes appearing at the quadratic order, 

• the space-time dependences are given by e^ ±k ( b > 1±k ( b ) 2 - l ' x ' l e _ ^ ±w< ' k ( b ) 1 - l± " ( ' k ( b ) 2 ^ t respectively. 

One can similarly work out that f( 2 ' 1,0 '(k, ki, fc 2 , x, x\) are dertermined by J 4.( 1 ) m ( Ile< i)' out ( n °q)(fc ; k^), the boundary 
conditions appearing at the linear order in the amplitude expansion. 



We can now specialize to the holographic Feynman propagator as given by (|10l| ). We can see that we reproduce the 
general form (|31) at the second order in the amplitude expansion. We get /( 2 '°>^"(fc, k±, /c 2 , x) similarly from the bound- 
ary conditions at the quadratic order and satisfy the properties listed above. Once again these satisfy the symmetry 
requirements but not reality constraints at all x and t unless both A ln ( neq )(k, fc(b)i, ^0)2) and A out ( neq )(&;, fc( b )i, fc( b ) 2 ) 
are set to zero, thus the solution is taken to be exactly the same as given by regularity at the horizon. Similarly 
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k2,x,xi) in (BlJ) should also vanish. This result generalizes to higher orders in the 



one can show that f^ 2 ' l,0 \k,ki. 
amplitude expansion. 

Thus we show that the holographic Feynman propagator is given by (105) to all orders in the ampl itude expansion. 
In other words only the solution which is regular at the horizon contributes in the linear response ( 101 ) holographically 
in determining the non-equilibrium Feynman propagator. 

Given that both the non-equilibrium retarded Green's function and Feynman propagator are determined by the 
regular solution, all non-equilibrium propagators are also determined as well by the regular solution in the gravitational 
response ( |98| ) . Thus the boundary conditions which determine the non-equilibrium corrections to the Green's functions 
are determined by regularity alone. 



D. On why regularity at the horizon is sufficient 



Our main result is that the regularity at the horizon and linear response theory together determine all the non- 
equilibrum Green's functions holographically. We can try to understand the physical significance of this result in 
terms of the dual field theory. 

Clearly the non-equilibrium dynamics of the expectation value of the operator, as for instance the energy-momentum 
tensor, is determined by regularity at the horizon. Indeed this is how hydrodynamic transport coefficients are deter- 
mined. In field theory this is given by the effective action T[0] as discussed in section II. 

We recall from section II. A that the non-equilibrium evolution of the propagators is given by the effective action 
T[0, G] evaluated over the Schwinger-Keldysh contour. In section II. A we also showed that the effective action T[0, G] 
can be obtain from T[0] applying functional identities. Thus, we should require no new information to determine the 
dynamics of non-equilibrium Green's functions other than those required to determine the non-equilibrium evolution 
of the expectation value of the operator. 

Holographically, this should imply that, just as regularity at the horizon is a sufficient principle to determine the 
non-equilibrium evolution of the expectation value of the operator, it should be a sufficient principle to determine 
the non-equilibrium Green's functions also. So the results of section II. A give a field theoretic justification of the 
holographic arguments advanced here. 



V. THE HOLOGRAPHIC NON-EQUILIBRIUM FLUCTUATION-DISSIPATION RELATION 



In the previous section we have proved that the holographic non-equilibrium Feynman propagator is given by (105). 
This result is true in non-equilibrium states close to thermal equilibrium, or more precisely in states perturbatively 
connected to equilibrium in derivative and amplitude expansions. This result holds in the absence of any external 
forces. Also this result is valid even when the non-linearities in the dynamics of hydrodynamic and non-hydrodynamic 
variables characterizing the non-equilibrium state are taken into account. 



From the holographic Feynman propagator (105), one can readily obtain the non-equilibrium statistical function 
using (|l8|). Thus we obtain the holographic non-equilibrium fluctuation-dissipation relation : 

G K {k,x) = i(2n BE (oj) + ljImG R {k,x) = -i(n BE {oj) + ^jA(k,x). (106) 

We recall again that the Bose-Einstein distribution 7t.be (w) above is determined by the temperature of final thermal 
equilibrium. Thus at long times we recover the usual equilibrium fluctuation-dissipation relation. 

We may wonder if this result is in contradiction with locality and causality of the dual field theory, because it 
depends on the temperature of thermal equilibrium to be attained in the future. This is not the case because of the 
following reasons : 



The fluctuation-dissipation relation (106) relates the Wigner transformed statistical function G/c(k, x) and 
A{k, x). The Wigner transform involves Fourier transform in the relative coordinate x\ — x-i both in Gjc(xi, x%) 
and A(x\, X2). The Fourier transform receives contributions from all values of the relative coordinate x\ — x%, 
thus both G/c(k,x) and A(k,x) are non-local in space and time by construction. 



• The result (106) is true strictly in the absence of external forces. In such situations, no energy is being pumped 
in externally, therefore the final equilibrium temperature is fixed by the total energy of the initial state. So 
there is no teleological adjustment involved as it may naively appear to be the case. 

Crucially G/c(k,x) and A(k,x) behave as local objects in space-time typically at weak coupling. In this case, the 
Kadanoff-Baym equations governing their evolution reduces approximately to Boltzmann equation which is local in 
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space and time. At strong coupling we do not expect this to be the case, therefore we have no reason to expect 
that Gjc(k,x) and A(k, x) will be related by local data. What we find is that the total conserved energy in the 
laboratory frame is sufficient to relate them - we do not need any more details of the non-equilibrium state, provided 
the field theory can be described holographically and the classical gravity approximation with a minimally coupled 
bulk scalar is valid. Individually Gx(fc,x) and A(k,x) though, carries detailed information of the relaxational modes 
of the system and their non-linear dynamics. 

We also note that in the late equilibrium state where our derivative and amplitude expansions are valid, the 
background hydrodynamic and non-hydrodynamic variables satisfy generic phenomenological equations which can be 
derived from classical gravity (for details see appendix B). These phenomenological equations, just like Navier-Stokcs 
equation hold irrespctively of initial conditions. The dual geometries have regular future horizons. Nevertheless, we 
can expect that not all solutions of these generic equations can be lifted to solutions in the field theory. Only a class 
of dual non-equilibrium geometries will also be regular in the past, and these will be genuine solutions in field theory. 
However, regularity in the past cannot be analyzed in perturbative derivative and amplitude expansions. Thus there 
will indeed be hidden constrants in the initial conditions for Su, 5T and 7rj" h ^ coming from the requirement of regularity 
in the past. Nevertheless, these do not affect our main conclusions, at least for the large class of non-equilibrium 
states which admit permit perturbative derivative and amplitude expansions in the future. 

We may also ask how this result may generalize to higher point non-equilibrium Green's functions. To address this, 
we have to study the backreaction of the scalar field on the background geometry with an arbitrary external source 
at the boundary. Regarding the generalization we make the following conjectures : 

• In the classical gravity approximation, the fully causal higher point Green's function @(t n — i„_i)....0(i2 — 
ti)([0(x n ), [0(x n -i), [0(x 2 ), 0(xi)] ]]j can be determined from the regular solution alone. 

• The other higher point Green's functions can be obtained from the causal higher-point function using their 
general field-theoretic properties and the consistency of their allowed forms with the classical gravity approxi- 
mation. 

The first point above simply says that a regular future horizon determines the full causal structure and hence the 
causal Green's functions. The second point above generalizes our derivation of the holographic Fcynman propagator. 
This conjecture says that we do not need to construct any partition function to obtain non-equilibrium correlation 
functions - general dynamical principles should suffice to determine them as in field theory. In future work, we would 
like to test or prove this conjecture. 



VI. DO OUR RESULTS SURVIVE STRINGY CORRECTIONS TO CLASSICAL GRAVITY? 

We will find here that we cannot certainly say that our results will survive the stringy corrections which typically 
give l/VA corrections in the strong coupling limit as a function of the dimensionless coupling A. In gauge theories, A 
is the 't Hooft coupling. 

Nevertheless all of our conclusions generalize even if there is stringy corrections in gravity dynamics (involving 
higher orders in curvatures and their derivatives), provided the scalar field dual to the operator is minimally coupled 
to gravity. The latter is certainly not guaranteed, but could be true in supersymmetric contexts. 

Firstly, we will show that our presciptions which give solutions which are regular at the future horizon in the 
non-equilibrium geometry, indeed remain valid as long as the bulk scalar field is minimally coupled to gravity. 

We observe that the structure of the homogeneous incoming solution near the horizon can be determined from 
geometric optics, if the bulk scalar field is minimally coupled to gravity. We can certainly construct an appropriate 
function of r which we denote as r* (r) such that the incoming radial null geodesic at the horizon is : 

v = t — r*(r). 

Clearly r„(r) has to increase monotonically as r moves towards the horizon because of blue-shifting. The homogeneous 
incoming wave solution in the thermal background by a quasi-normal mode perturbation will always behave like : 

py e ~ i ("+"(b)) t ' 

near the horizon at the leading order, as the geometrical optics approximation is always good at the horizon due to 
the blue- shifting. Therefore, as long as thermal geometry is stable, the integration over co^) which we need to do to 
put the metric on-shell, will produce a divergent factor : 



e -Imw (b) (k (b) )r*(r)_ 
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If the thermal background is stable, all it's quasi-normal fluctuatons should satisfy Imw^jt^jb)) < 0. Therefore the 
homogeneous incoming solution will always be divergent at the horizon. 

On the other hand the argument that the outgoing modes will cause divergent backreaction at the horizon is based 
on analyticity in u - this remains true as long as the thermal backgound is stable. 

In any consistent classical gravity, the black brane dual to the thermal background should be stable, thus we need 
to cut down both the homogeneous incoming and outgoing modes for regularity at the horizon. 

The chain of arguments determining non-equil ibriu m propagators, similarly will also carry through all the way to the 
non-equilibrium fluctuation-dissipation relation ( |106| ). This will be true provided the black brane do not support quasi- 
normal modes where both the momentum and frequency could be imaginary, or where the frequency/momentum is 
imaginary while the momentum/frequency vanishes respectively. One can check, if this happens, the reality arguments 
which cut down other boundary conditions to determine non-equilibirum propagators other than those required by 
regularity, will no longer hold. 

In fact, a quasinormal mode with imaginary momentum at zero frequency, will signify the existence of hair. Also, a 
quasinormal mode with imaginary frequency at zero momentum can be interpreted as a thermal tensor-like Goldstone 
mode. The first situation is unlikely as usually black branes do not have hair. The second situation is unlikely as 
usally in field theories we do not have spontaneous breaking of Lorentz symmetry. 

Even if the bulk scalar remains minimally coupled taking stringy corrections into account, we still need some 
care. Though the equations of motion may admit systematic perturbative expansion in 1/V% all solutions which 
thermalizc to black branes need not admit such expansion. Also we cannot take the limit A — > smoothly, as in this 
limit curvature diverges asymptotically. Therefore, we still need to be on the strong coupling side, for our results to 
generalize. 



We may conclude that the nonequilibrium fluctuation-dissipation relation (106) requires not only the validity of the 
classical gravity approximation, but also the classical theory of gravity to be nearly Einstein's theory, and the bulk 
scalar field to be minimally coupled to gravity. 

VII. CONCLUDING REMARKS 

The key points of this paper may be summarized as follows. 

• The holographic non-equilibrium spectral function carries systematic information about the relaxation modes 
of the system and their collective non-linear dynamics. 

• The holographic prescription for obtaining the non-equilibrium spectral function is independent of the non- 
equilibrium state i.e. the gravitational background and is obtained from requiring regularity at the horizon at 
strong coupling and large N. 

• This holographic prescription should also hold in any theory of classical gravity which has a stable thermal 
background as a solution as long as the bulk scalar dual to the bosonic operator is minimally coupled to gravity. 

• It is possible to map the parametrization of non-equilibrium Green's functions required for their consistent 
derivative/amplitude expansions to linear response functions in gravity using arbitrary boundary conditions for 
the non-equilibrium modes. 

• Field-theoretic consistency determines the holographic non-equilibrium Feynman propagator. 

• From the holographic non-equilibrium Feynman propagator, we obtain the non-equilibrium fluctuation- 
dissipation relation. 

• The holographic non-equilibrium fluctuation-dissipation relation holds universally for any non-equilibrium state 
at strong coupling in the classical gravity approximation and thus provides a potentially strong test of the 
applicability of holographic duality. 

Except for the first two points, all these have been derived in this paper. 

We may now discuss how can the holographic scenario be tested in RHIC or ALICE. Our discussion will assume 
that the plasma produced by the heavy-ion collision undergo three stages of evolution - (i) the initial phase (ii) a 
strongly coupled nearly conformal phase, and (iii) the final hadron gas phase. The classical gravity approximation 
studied here can only apply to the middle phase of evolution. This picture is supported by lattice studies as in |}l[] 
where it was found that QCD is inded close to a strongly coupled fixed point at a scale of about 175 MeV (which is 
nearly the temperature of the quark-gluon plasma produced at RHIC/ ALICE) for small chemical potentials. This is 



also consistent with the observed incredibly short time of thermalization (« 1 fm/c) 32, p3|. 
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We note that in the early stage of the ultra-relativistic collisions we may use perturbative QCD which gives a 
quantum kinetic theory involving quarks and gluons. In this case, we may use models like |j34j| . We may also use color 
glass condensate models where the two particle irreducible (2PI) effective action technique has been used to obtain 
growth of fluctuations and correlations of various operators [|3^] . In the late stage of expansion of the fireball, we can 
use the Ultra-relativistic Quantum Molecular Dynamics (UrQMD) simulations which models the quantum kinetics of 
the hadron gas [p6| . The holographic fluctuation and correlation of the chiral condensate needs to be matched with 
the quantum kinetics of the early and late phase of the expansion to obtain a complete picture that can be tested 
experimentally. 

In the holographic phase, we may use the gravitational backgrounds provided by Einstein's gravity with anti de-Sitter 
boundary conditions to calculate the non-equilibrium spectral function. The non-equilibrium spectral function gives us 
a way to determine the collective non-linear dynamics of the relaxation modes and thus check if the phenomenological 
equations obtained from gravity really describe the space-time evolution of the fireball in the strongly coupled phase. 
The latter can also be tested independently by detailed measurement of transverse momentum spectra and elliptic 
flow coefficient of the hadron gas. A mutual confirmation will be a strong test for holography. 

The non-equilibrium spectral function combined with the non-equilibrium fluctuation-dissipation relation gives 
sufficient data to match with correlations of emission of pions and their resonances in the later stage of expansion, 
as pions and their resonances have the same quantum numbers as the chiral condensate. Thus we can build a more 
sophisticated theory of pion interferometry to reconstruct the expansion of the fireball and test the holographic 
non-equilibrium fluctuation-dissipation relation. As the fluctuation-dissipation relation holds universally for any non- 
equilibrium state in the classical gravity approximation, it's validation will indeed be a very strong test for holography. 

In fact, it has been observed that statistical models using pure thermal estimates work very well for predicting 
the spectrum of hadrons produced from the fireball with the temperature being set to the value at which the fireball 
thermalizes fl37| , [[}2| . This is surprising given that the phase transition does not occur adiabatically. This may already 
hint that the holographic non-equilibrium fluctuation dissipation relation is valid to a good accuracy. 

The statistical function of the chiral condensate carries direct information about the production of pions and their 
resonances. Thus a time-resolved study of emission of pions and their resonances will allow us to determine the 
non-equilibrium statistical function once we know how to match it's evolution to the early and final phases. 

A more refined time-resolved study of correlations of emitted pions and their resonances will allow determination 
of the non-equilibrium spectral function. This also needs to be matched with the early and late stages of evolution. 
Having thus determined the non-equilibrium statistical and spectral functions individually, we can check the non- 
equilibrium fluctuation-dissipation relation. 

The matching can be done by noting the following that in the holographic phase, the non-equilibrium state, and 
also the non-equilibrium spectral and statistical functions, can be parametrized by a few non-equilibrium variables 
like the hydryodynamic variables and the shear-stress tensor. These variables represent the expectation value of a 
handful of operators like the energy-momentum tensor. Thus the matching can be realized by tracking the evolution 
of these operators. Some interesting progress has already been made in this direction Q . 

It thus looks possible that the holographic non-equilibrium fluctuation-dissipation relation can be tested in the 
future in heavy-ion collisions. This may also lead to an accurate experimental determination of the temperature at 
which the plasma thermalizes. 

It will be also interesting to understand how the non-equilibrium fluctuation-dissipation relation be measured in 
order parameter correlations in quantum critical systems. 
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Appendix A: Useful identities 

Consequences of symmetry and transposition 

Consider a Hermitean scalar operator O(x). The retarded Green's function for this operator is defined as 

g r { Xi ,x 2 ) = -ie(t 1 -t 2 ){[o(xi),o(x 2 )]y 

The advanced Green's function is defined as 

G A {xi,x 2 ) = iQ(t 2 - [0( Xl ), 0{x 2 )\ ). 



It follows from the definitions that 



Ga(xi,x 2 ) = Gr(x 2 ,x\). 



We may also see that 



G R * (xi,x 2 ) = G R (xi,x 2 ), G a * (xi,x 2 ) = G A (xi,x 2 ). 



(Al) 



(A2) 



(A3) 



(A4) 



Let us define the center-of-mass coordinate as x — (x± +x 2 )/2 and the relative coordinate r = x\ — x 2 . The Wigner 
transformed retarded and advanced Green's function are : 



Gn, A (k,x) = j d 4 re ,tT G M ( S ,r). 
It follows from ( |A^) that Ga{x,t) = Gr(x, —r). We thus see that 

G R (-k,x) = J d A re- lk - r G B {x,r) 
= [ d 4 re lk - r G R (x,-r) 



(A5) 



d 4 re tk - r G A {x,r) 
= G A (k,x). 

In the second lin e above we have changed variables from r to — r. 

Similarly, (A4) implies that Gr * (x,r) — Gr{x^t) and Ga * (x,r) = GA(x,r). It follows then 

G R *(k,x) = I d A re- lk r G R *{x,r) 

d^re-^GRix.r) 
G A (k,x). 



(A6) 



(A7) 



In the secon d li ne above, we have used Gr * (x, r) — Gr(x, r). Observing that the second line coincides with the first 
equality in (A6), our final conclusion follows. We can similarly prove that Ga * (k,x) = GR(k,x). 

We now turn to the Feynman propagator Gf(xi,x 2 ) which is symmetric. Therefore, Gf(x, r) = Gf(x, —r). 
Therefore it follows that the Wigner transform of the Feynman propagator should satisfy : 

G F {-k,x) = [ d A re~ lk - r G F {x,r) 

d 4 re lk - r G F {x,-r) 

d\e tk - r G F (x,r) 

(A8) 



G F (k, x). 
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In the second line above, we have changed variables from r to — r, and in the third line we have used Gf{x,t) 
G F (x,-r). 

Similar ly it follows from symmetry in x\ and X2 prior to Wigner transform that G/c(k,x) = G;c(—k,x). 



Using <\A§, we also note that ReG R (k,x) = (l/2)(G fi (fc,x) + G A (k,x)) = (l/2)(G A (-k,x) + G R (-k,x)) = 



Re Gr{— k, x) 
The statistical function is purely imaginary 

According to (|^), the inverse Wigner transform of the statistical function is given by the anti-commutator as 
below : 

G K (x,r) = -i({0(x + l),0(x-l)}). (A9) 



Clearly, 



Therefore, 



G K {x,r) = G K (x,-r) ) Gk, * (a;, r) = —Gjc(x, r). (A10) 

G K *(k,x) = / d A re- ik - r G^*(x,r) 
d 4 re~ zk - r G K {x,r) 



d 4 re* k - r G K ( Xl -r) 
= -J d 4 re lk - r G K (x,r) 

= -G K (k,x) (All) 

In the second line above, we have used Gk. * {x, r) — —Gk(x, r), in the third line we have changed variables from r to 
— r, and in the fourth line, we have used Gk(x, r) — Gk(x, —r). 

Clearly, Gk * (k, x) = —Gjc(k,x) implies that Gic(k, x) is purely imaginary. 

The spectral function and the retarded Green's function 

According to (|l|), the inverse Wigner transform of the spectral function is given by the commutator as below 

A(x,r) = ([o(x+ r -),o(x- r -)]). (A12) 

It follows from ([Al]) and(|Aj|) that 

A(x > r)=i(G R {x,r)-G A (x,r)^ = -2lxaG R {x, r), (A13) 
using 6(ti — £2) + 9(*2 — t%) = 1- Clearly then after Wigner transform, 

A(k,x) = i(G R (k,x)-G A (k,x)\ = -21mG R (k,x). (A14) 
Using (|Aq), we readily see from the above that 

A(k,x) = -A(-k,x). (A15) 



The Feynman propagator as a sum of the statistical and spectral functions 
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The definition of the Feynman propagator is : 

G F (x 1 ,x 2 ) =-i(T(0(xi)0(x 2 ))), (A16) 

where T denotes time ordering. 

It follows from the definition of the Feynman propagator that 

G F (x 1 ,x 2 ) = -i(0(xi)0(x 2 )J if h > t 2 

= ~(0(x 1 )0(x 2 ) + 0{x 2 )0[xi)) - l -(0{x 1 )0{x 2 ) - 0{x 2 )0{x x )) if h > t 2 

= G/c(xi,x 2 )-~A(x 1 ,x 2 ) if ti >* 2 . (A17) 

Similarly, 

G F {x 1 ,x 2 ) = -i(0{x 2 )0{x 1 yj if < 2 >h 

= ~(0{xi)0{x 2 ) + 0{x 2 )0{xi)) + l -(o[x x )0{x 2 ) - 0{x 2 )0{xxj) if t 2 > h 

= G lc (x 1 ,x 2 ) + ^A(x u x 2 )i£t 2 >t 1 . (A18) 



Combining these, we obtain 



2 

G F (xi,x 2 ) = G)c(xi,x 2 ) - -A(x-l,x 2 ) sign(ii - i 2 ). (A19) 



Appendix B: The general phenomenological equations for irreversible processes 

The phenomenological equations generalizing hydrodynamics which can be obtained from gravity can be argued to 
be as follows. We can do a kinematic Landau-Lifshitz decomposition of the energy-momentum tensor as below : 

V = eu v u v + P p n» + V«" ( B1 ) 

with u M being a time-like vector of unit norm (so u^u u = —1), and P^ v being the projection tensor in the spatial 
plane orthogonal to u 11 given by P^v = Up,u v + ?/ M „ . 

We can define a local temperature using the equation of state e(T) locally. We can also set p by using the local 
equation of state p(T). Then 7r MI , denotes the local non-equilirbium part of the energy-momentum tensor. 

Also is defined to be the local velocity of energy transport, therefore 

= 0. (B2) 

Thus 7r M „ has six independent components only. Together with u M and T, it gives ten independent variables to 
parametrize ten independent components of t^ v . 
Furthermore, conformal invariance requires 

e = 3p, ir^rf" = 0. (B3) 

This reduces the number of variables to nine, cutting down the independent components of tt^ from six to five. 

To proceed further, we need to decompose 7r Ml/ into two parts, a purely hydrodynamic part tt^J and a non- 
hydordynamic part 7rjiy which is dynamically independent of the hydrodynamic variables though they do couple to 



them M, [El . Thus, 



_ (h) , (nh) 



(B4) 



The fluid/gravity correspondence tells us that at strong coupling, purely hydrodynamic states exist in the dual field 
theory in the derivative expansion. Thus we can consistently set the non-hydrodynamic part 7r^" h ^ to zero. 
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The general form of the purely hydrodynamic part tt^J can be readily obtained order by order in the derivative 
expansion by constructing the conformally covariant tensors which are algebraic functionals of the derivatives of the 
hydrodynamic variables. The derivative expansion parameter is typical length scale of variation with respect to the 
thermal wavelength. The transport coefficients are determined via regularity at the future horizon p8fl. Thus 



rOO 



T)P/P v a (d p u a + d a u p - ^P pa {d ■ u)) + 0(e 2 ). 



(B5) 



Putting 7r^ nh ' > 



in the form of t pv in (Bl), and demanding d^t pu = results in the hydrody- 



0, plugging - 
namic equations of motion. 

Energy and momentum are always conserved and this follows directly from the constraints of Einstein's equations. 
This implies that 



d» V - d» (e(T) u p u v + p(T) P pv + tt£) + ^) = 0. 



(B6) 



The above are sufficient to determine the evolution of w M and T, but not that of the non-hydrodynamic shear-stress 
tensor ir^K We require five more equations to determine the latter. These should be provided by the regularity of the 
horizon in gravity as in the case of homogeneous relaxation. We need a separate expansion parameter - the amplitude 
parameter /p. As in the case of homogeneous relaxation, though the spatial derivates of ^if^ can be expected 
to be small near equilibrium, the time-derivatives will be 0(T). The Lorentz and Weyl covariant generalization of 
the time-derivative is T> = (u ■ d) + .... 

The equation of motion for Trj^ should be such that (i) it can be consistently put to zero (ii) it is Lorentz and 
Weyl covariant and (iii) all time-derivatives are summed over at each order in the amplitude expansion. Requiring 
these, it follows that this equation should take the form EM, p5[: 



\n=0 I 



n— m— 

+0(e 2 5, eS 2 ,S 3 ) 



(B7) 



up to given orders in the derivative/amplitude expansion. Above (...) and [...] denotes symmetrization and anti- 
symmetrization of enclosed indices respectively. 

Let us see how we can derive these equations from gravity for the special case of homogeneous relaxation. In this 
case, the flow is at rest, so we can go to an inertial frame where u M = (1, 0, 0, 0). Also the temperature T is constant in 
space and time. Furthermore (B2) implies that the non-zero components of 7r M „ are the spatial components TTij. The 
conservation of energy and momentum implies that TTij(t) is an arbitrary function of time. Nevertheless, regularity at 
the horizon can be expected to be guaranteed only when mj (t) follows a definite equation of motion. 

One can construct the dual metric in such cases perturbatively in the so-called amplitude expansion 7Tj-" h '(i)/p. 
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Omitting the details, the metric up to second order in the amplitude expansion, takes the form below pi| : 
I 2 dr 2 I 2 



ds z = 



/(=?) 

OO j 



71=1 

oo n 



+ E E ^ fi M (£) (|) ro 4f ) 

n— 1 m— 1 

oo n j j 

E E ^ + " + "7f - m) (£) (|R h) (|) m 4f ) ^ 



n— 1 m—1 

oo n 



'» + " + '"^" (i) --"(i) + <• «■» - Nl) 4f» (|R» 

n— 1 m—1 \ / 

+0(<5 3 )). (B8) 

We note that at each order in the amplitude expansion we have summed over all time-derivatives of 7r£° (t) ■ Thus at 
the first order in amplitude expansion, we get an infinite number of radial functions /(b«) correponding to the n-th 
time-derivative. These can be determined uniquely by requiring /^'^(s) — s 4 + 0(s 8 ) and /^'"•'(s) = 0(s 4+Tl ) for 
n > 1. In these Schwarzchild coordinates it turns out that /t 1 '™) vanishes for odd n. Similarly at the second order in 
the amplitude expansion we get infinite number of radial functions f( 2 ' n ' m ^ j-( 2 - n > m ) anc j j^ n , m ) a vj j wn j c j 1 can \y e 
determined uniquely from their boundary behavior f^ 2 ' n -' m ^ = 0(s 8+n+m ) for i = 1,2,3. Also in these Schwarzchild 
coordinates j| 2 ' n,m ) vanishes when n + m is odd. 



The regularity analysis of the metric ( |B8| ) is subtle and involves it's translation to Eddington-Finkclstcin coordinates. 
One finds that the metric is regular, if up to second order in the amplitude expansion tt^J (t) satisfies the following 
equation of motion [El] : 



n—0 7i—0 m—0 

+0(6 3 ) = 0. (B9) 

In order to see regularity of the metric at future horizon, order by order in the amplitude expansion, it is necessary to 
sum over all time-derivatives at each order. One can determine the non-dydrodynamic phenomenological coefficients 
£)(!,«) anc [ jj(2,n,m) j n t erms G f complicated recursion relations. The first few terms are : 

= —ttT, D ™ =-(|4 lu2) (7 rT) 2 , etc, = —L^, etc. (BIO) 

That it is necessary to sum over all time-derivatives at each order in the amplitude expansion to obtain regularity at 
the horizon, can be understood from the first order in the amplitude expansion itself. It can be shown that the series 
Dft' 71 ^ (~iuj) n has simple zeroes at the location of the discrete homogenous quasi-normal modes (with k = 0) in the 
complex lower half plane. The location of such quasi-normal modes are approximately given by [ p9[ 

w (n) «7rT ± 1.2139- 0.775 i ± 2n(l T i) ■ (Bll) 

We observe that both the real and imaginary parts are of 0(T) signifying that the time-derivatives of 7r[" h ' 1 are of 
0(T). Thus, we do need to sum over all time-derivatives as we cannot do a small frequency expansion even in the 
linearized approximation. 



The equation (B7) reduces to flB9j) in the special case when T and u M are constant is space and time, and furthermore 



we go to the laboratory frame where u M = (1,0,0,0). Thus D^' n ^ and D^' n,m > should be as given in (B1C). The 



coefficients X^ and X^ denote coupling of ir^J 1 ^ and it's local time-derivatives to hydrodynamic variables. To 
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determine these, one needs to construct metrics with regular future horizons for configuration where and T vary 
spatially and temporally, while 7r^ I " h ' 1 is non-zero. 



Appendix C: Non-equilibrium shifts in quasiparticle dispersion relations 

Holographically the vanishing of the source gives the quasi-particle poles (which can be so broad that we may not 
call them quasi-particles) . At equilibrium J m ^ cq -\uj, k) vanishes only when u> takes certain discrete values at a given 
k. We may select a branch of quasi-particle pole given by: 

= ,,,( c( lh\c) surt, that ^in(eq)/ (eq) 



W = W 



(k) such that J" n ^(w^(k),k) = 0. (CI) 



These correspond to the poles which can be complex generally. The non-equilibrium modification depends on space 
and time, and takes the form 

w = w(° q )(k) + fe(k,x,t). (C2) 

The non-equilibrium shift in the pole <5w(k, x, t) can be obtained by solving J7(x, t) — perturbatively |T^ |. This 
amounts to solving the linear complex equation : 

fe(k,x,t)9 wl 7 in ( oq) (w = u; (cq) (k),k) = -J in ( ne ^L=w (eq) (k),^ 

-^W L = u W(k),k, -k (b] je <BflU w (k m ))i e I,n ''w (k w )( e-*( b ) ,3C . (C3) 

The non-equilibrium parts of the source, or equivalently the right hand side of the above is determined completely by 
the equilibrium source by our prescription at the horizon. This is precisely what we expect in the dual field theory. 
The non-equilibrium shift in the pole should be determined solely by the infrared behavior. For example, if the 
temperatur e flu ctuates, so should the thermal mass. This gives a good justification of our prescrition. Our general 
expression (C3) shows that the pole of the resummed non-equilibrium propagator can depend also on the velocity and 
shear-stress perturb atio ns also. 

It is evident from ( |C3| ) that 8uj (k, x, t) will be comoving with the relaxational mode constituting the non-equilibrium 
state (or the dual quasinormal mode) as below : 

<Jw(k,x,i) = < 5o;(k,k (b) )e- lRo "( b ) (k ( b ) ) *e Im "( b ) (k ( b ) ) *e' ik ( b )- x 

+Soj{k, _k (b) )e 4RoLJ ( b ) (k ( b ) ))t e Im "( b > (k ( b ) )t e- 4k < b )- x . (C4) 

Furthermore, <5w(k, kn^) will have a consistent expansion in the derivative/amplitude expansion which follows from 
that of J in ^((J,k,k(b)). For example, in the case of the hydrodynamic shear wave background 6w(k, k(h)) takes 
the form : 

(5w(k, k(h)) — Scoi(k)5u ■ k + Su2(k)kikjk^iSuj (C5) 

up to first order in the derivative expansion. We thus explicitly see in the above example how our holographic 
prescription gives the non-equilibrium shift in the quasi-particle pole parametrized in terms of the background velocity 
perturbation. 
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Appendix D: Wigner transform of the holographic retarded Green's function 

Writing down explicitly we get, 
G Jl (x 1 ,f 1 ,x a ,f 3 ) = J dud 3 k^^lA = J du J d z k (Dl) 

I 0( c i)( tj j i k)e~ iwtl e ik ' Xl + O in ^ noq ^(a;, k, k( b ))e _ ^" +Rc " (b ^ k(b '^* 1 e ImLJ(b '^ k(b ^ tl e^ k+k(b) ' ) ' Xl 



_^Qin(aeq) ^ _j < .^^ e -i(w-Rc w (b) (k (b) ))ti glm w (b) (k (b) )ti e i(k-k (b) )-xi ^ j ^ J r (°q)( W; k) 
_|_ j-in(noq) ^ ^ ^ g -i(a;+Re w (b) (k (b) ))t 2 glm w (b) (k (b) )t 2 gi(k+k (b) )-x 2 
_|_ j-in(ncq) ^ _ ^ ^^ g -i(a;-Re u (b) (k (b) ))t 2 glm w (b) (k (b) )t 2 gi(k-k (b) ) -x 2 ^ _ 

The above can be approximated in the derivative/amplitude expansion as below : 

( O in ( ncq Hw,k,k (b) ) ., . n v t 

1 i v 1 ' Wj, ptk(bV x i r -»Rc"(h)(kfh))«i r Imcj lb )(k lb i)ti 

^ + 0(«*i>(w,k) 

/nin(neq)/, , i_ i_ \ 

^ *(b); -tk (b) -xi flUa; (b) (k (b) )ti Imh> (b) (k (b) )ti 

+ eX^.k) 



g-«*>t 2 g ik -x 2 



,7ta(neq)( W)k>k ) 



( b )^ c iU. lb ) x 2 g-iRo w (b) (k( b) )t 2 glm w (b) (k (b) )t 2 

J"(«0(w,k) 

/7-in(neq)/, . l_ l_ \ \ 

O \"j*2 ^(b); e -ik (br x 2e iRec t;(b) (k fb) )t 2e Ima; (b) (k (b) )t 2 ( D2 ) 

j( cq )(w,k) ' r 



It is easy to take the Wigner transform of the above to obtain (B4 
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